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CHAPTER I 
INTRODUCTION 


Flush mounted antennas for the space shuttle vehicle will 
generally consist of dielectric covered slots in a metallic surface 
where the dielectric cover is employed as a protective heat shield 
for the re-usable space vehicle. A typical local antenna geometry 
is illustrated in Fig. 1 where the truncated dielectric cover is 
shown recessed and flush mounted in a perfectly-conducting surface. 

The work described in this report deals with the analysis of the 
electromagnetic radiation from such an antenna configuration. In 
particular, it is the purpose of this analysis to assess the effects 
of the structural discontinuity caused by the truncation of the 
dielectric cover (at A and B in Fig. 1) with a view towards pro- 
viding analytical and numerical results which may be useful in the 
design of such antennas. 

In the present work, the configuration illustrated in Fig. 1 is 
two dimensional, and the infinitely long axial slot excites only TM 
type electromagnetic fields; here TM implies that the magnetic field 
intensity is entirely parallel to the axial slot (i.e., it is directed 
normal to the plane of the paper). This configuration supports 
surface waves excited by the slot (in addition to the fields directly 
radiated by the slot). In the present analysis, the thickness, b 
and the relative permitivity, e r of the dielectric cover are chosen so 
that only the dominant TMq surface wave is excited for a given op- 
ating frequency. The radiation pattern of this antenna can be con- 
veniently calculated via the geometrical theory of diffraction[l ]; 
however, the diffraction and reflection coefficients associated with 
the truncation of the dielectric cover must be found. Expressions 
for these coefficients are derived in this report. 

In terms of the geometrical theory of diffraction, the total 
far-zone field is composed of a geometrical optics field, and a dif- 
fracted field. The geometrical optics and diffracted fields are 
indicated in Fig. 1. The geometrical optics field for this con- 
figuration is the direct radiation from the slot to the field point P 
along the ray path s 0 ; whereas, the surface waves excited by the 
slot impinge on the terminations at A and B to produce singly dif- 
fracted waves which arrive at P via the ray paths s-| and S 2 * The 
surface waves excited by the slot also produce reflected waves at 
the terminations. These reflected surface waves are in turn inci- 
dent on the opposite terminations thereby producing additional re- 
flected waves, and also diffracted waves which arrive at P along 
the ray paths s-j and S2» and so on. The diffracted field at P 
is therefore the sum of the fields associated with the singly 
diffracted and multiply reflected-diffracted rays arriving at P. 

The latter contribution can be summed into a closed form using a 



self-consistent method. The geometrical optics field and the field of 
the surface waves (incident on the terminations) are obtained from the 
solution to the problem of a narrow slot radiating through a dielectric 
covered ground plane of infinite extent. .The diffracted and reflected 
fields produced by the terminations at A and B are described in terms 
of the appropriate surface wave diffraction and reflection coefficients 
just mentioned. 



Fig. 1— A typical antenna configuration involving a dielectric 
covered slot in a perfectly conducting surface. 1 


It is apparent that the superposition principle can be used to 
extend this analysis for the single narrow slot directly to an 
array of slots radiating in the presence of a dielectric cover, 
or to an extended aperture in the presence of a dielectric cover.'' 

Figure 2 illustrates the canonical problem pertinent to the ; r 
analysis of the terminated dielectric cover shown in Fig. 1. The 
semi-infinite dielectric slab recessed in a perfectly-conducting 
surface is excited by a TM 0 surface wave which is incident from ‘ 
the left on the dielectric termination at z = 0. In this report, 
three methods are employed to solve the above canonical problem; 
they are: 

T. a formally exact solution using the generalized scattering 
matrix technique combined with function theoretic methods, 
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2. an approximate solution using the variational technique, 

,'and - ;• . - 

0 ; * 

3.. an exact solution of an approximate model wherein the 

boundary conditions pertaining to the grounded dielectric slab 
are replaced by an equivalent surface reactance for x = 0 and 
z < 0 as illustrated in Fig. 3. One notes that this model also 
approximately describes the diffraction of a surface wave on the 
complementary structure of Fig. 4. It is conjectured that the ter- 
minated surface impedance model (Fig. 3) is a better approximation 
for the dielectric slab on the ground plane (Fig. 4) than for the 
recessed dielectric slab configuration of Fig. 2, because the latter 
structure has an additional discontinuity in the form of the step in 
the ground plane. This is evident by noting that as the relative 
permittivity of the dielectric slab approaches unity, the 
former structure reduces to the smooth ground plane, whereas a 
residual step is present in the latter structure. In Chapter IV 
it is shown that this conducting step significantly affects the 
reflection coefficient quite apart from the terminated dielectric 
slab. 


The. boundary value problems associated with the configurations 
shown in Figs. 2 and 3 are not ammenable to solutions by the 
separation of variables technique due to the complicated nature 
of the boundary conditions. The configuration in Fig. 3 is identi- 
fiable as the junction of two semi-infinite sub regions thereby 
suggesting a solution based on the Wiener-Hopf type factorization 
procedure. Thus, the two-part boundary value problem of Fig. 3 is 
formulated in terms of dual integral equations that can be solved 
exactly using the method of factorization. The factorization for 
this problem is achieved via an appropriate limiting operation on 
the Wiener-Hopf factors associated with a related closed region 
boundary value problem (a closed region problem is one for which 
no radiation takes place); this technique was originally used by 
Bates and Mittra[2] in their Wiener-Hopf solution, to the problem 
of waveguide excitation of dielectric and plasma slabs. Such a 
factorization procedure allows one to obtain factorization functions 
which are generally far more convenient, for numerical processing 
than those obtained by a formal factorization scheme. The details 
of this solution are presented in Chapter II. This solutipn yields 
the expressions for the surface wave reflection and diffraction 
coefficients for the approximate configuration in Fig. 3. Several 
authors have treated two-part boundary value problems whose solutions 
contain the solution to this problem as a special case; however, 
most of these authors do not provide numerical results, and where 
numerical results are given, they are only for the magnitudes of. 
the quantities of interest. One of the reasons for this is the 
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X 




( X = 0,Z<0 ) ( X = 0 , Z > 0 ) 

Fig. 3 —An approximation of the configurations 
shown in Figs. 2 and 4. 



DIELECTRIC SLAB PLANE 


Fig. 4 —A configuration related to that in Fig. 2. 
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complicated nature of their solutions. Our primary aim in treating 
the two-part boundary value problem of Fig. 3 is to obtain solutions 
'■ 'to quantities of interest in a form which is convenient for cal- 
culation and subsequent application to antenna design. Only one 
reference (Weinstein [5]) appears to provide analytical results which 
are tractable for calculation', in this case, the values obtainable 
from the present results and from those in [5] (which are based on 
a different factorization scheme) provide a valuable mutual check. 

There appears to be an error in the expressions for the surface wave 
reflection and diffraction coefficients given in [5]; if one changes 
the final results in Reference [5] to make them consistent with the 
development which preceeds the final results, it is found that 
the numerical results calculated via [5] agree well with the calcu- 
lations based on the results derived in this report. In particular, 
the values for the magnitudes of the reflection and diffraction coeffi- 
cients calculated by these two methods agree perfectly, whereas, there 
is only a slight difference in the corresponding values for the phase 
of these coefficients. 

The complexity of the solutions and expressions for the parameters 
of interest which are given by Kay [3] and Maliuzhinets[4] appear to be 
fairly representative of what one might encounter when consulting refer- 
ences dealing with two part boundary value problems; their expressions 
for the phase of the reflection and diffraction coefficients contain 
integrals which appear intractable for numerical computations; in fact, 
the results for both the amplitude and the phase of these coefficients 
. available in Reference [4] appear intractable for numerical computations 
(when the results are specialized to our case). 

Angulo[6] has given an approximate solution to the problem il- 
lustrated in Fig. 4. Using the continuity of the tangential com- 
ponent of the magnetic field in the aperture z = 0, x >_ 0, Angulo ob- 
tained an integral equation for the equivalent magnetic current. He 
employed the variational method to obtain a stationary expression for 
the aperture impedance of the incident surface wave. However, in his 
solution, he approximated the aperture distribution (magnetic current) 
with only the surface wave field. In related problems treated here we 
show that such a trial function is too crude to obtain accurate values 
of the surface wave reflection coefficient. It is clear that the radi- 
ation pattern can be calculated from this single term trial function 
without using the variational method; these patterns are found to be 
reasonably accurate, except at the lower power levels. 

The configuration in Fig. 2 has a transverse discontinuity (due 
to the step) in addition to the longitudinal discontinuity which 
defines the junction of two semi-infinite sub regions (in the longi- 
tudinal direction); hence, this problem gives rise to a modified 
Wiener-Hopf equation. Instead of solving the modified Wiener-Hopf 
problem, the solution to the problem in Fig. 2 is built up from 
the solutions to two simpler problems which can be treated via 
the Wiener-Hopf technique. Specifically, the solution to this 
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problem is formulated using a generalized scattering matrix[7], 
in which the elements of the generalized scattering matrix are 
deduced from the Wiener-Hopf solutions to two appropriate auxiliary 
canonical problems. Both auxiliary problems give rise to identical 
Wiener-Hopf factorization functions. It was indicated previously * 
in the discussion dealing with the solution to the problem in 
Fig. 3 that the- Wiener-Hopf factors for the open region could be 
found from a limiting operation on the Wiener-Hopf factors for the : 
related closed region, and that this is to be generally preferred . / 
over the formal factorization procedure for reasons of convenience 
in calculating the factors; unfortunately, the Wiener-Hopf fac- 
torization function associated with these auxiliary canonical ■ 
problems do not yield a convergent representation for the factors • e 
obtained via the limiting procedure (see Appendix I, Section III). 

It is therefore necessary to resort to a formal scheme in order to t' 
factorize the Wiener-Hopf Kernel for the auxiliary canonical jv." 

problems (for details, see Appendix II). It should be noted the 
generalized scattering matrix formulation leads to a formally , 
exact solution for the boundary value problem corresponding to 
Fig. 2. The results for the reflection and diffraction coefficients 
obtained from this solution are expressed in terms of appropriate 
scattering matrices which are of infinite order; however, accurate 
computations are possible by truncating the size of the matrices so 
that one deals with finite order matrices. For a given dielectric 
constant, operating frequency, and slab thickness, the size of the 
finite order matrices is determined by the accuracy desired. 

Numerical calculations indicate that in most practical cases, 
sufficiently accurate results may be obtained by dealing with a 
relatively small matrix. A discussion of the auxiliary canonical 
problems, their solutions, and the manner in which these are used 
to determine the generalized scattering coefficients appearing in 
the expressions for reflection and diffraction coefficients, is 
presented in Chapter III. In addition, an approximate solution to 
this canonical problem which is based on the variational technique 
is given in Chapter III. 

The use of the surface wave diffraction and reflection coeffi- 
cients to analyze the radiation from the antenna configuration of 
Fig. 1 via the GTD is described in Chapter IV. Chapter IV also 
contains numerical results for the reflection and diffraction 
coefficients which are derived in Chapters II and III. By a 
comparison of the numerical results, it is concluded that the 
approximate physical model in Fig. 3 is not adequate for analyzing 
the configuration in Fig. 2 unless the conducting step in Fig. 2 
is extremely small (much less than 1/50 x^, where Xj = the wave- 
length in the dielectric); hence, the formally exact solution 
given in Chapter III is to be preferred even for relatively thin 
dielectric slabs. The technological relevance of this research to 
the design of arrays radiating in the presence of a dielectric cover 
and to the design of surface wave antennas is pointed out in 
Chapter IV. 
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An e~ 1ut time convention is used in this report; this time 
convention is generally employed. by physicists; whereas, the. 
electrical engineers prefer an e la) * time dependence. The e”i“t 
time dependence is employed here because it is commonly used for 
solving problems by the Wiener-Hopf method, and the notation and 
definitions employed are directly connected with this time con- 
vention. 


The authors have benefitted from the helpful discussions with 
Professors H. D. Colson, J. H. Richmond and L. Peters, Jr. Special 
thanks are due to Professor D. C. Chang of the University of Colorado 
for pointing out the utility of the generalized scattering matrix 
technique to the authors. 
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CHAPTER II 

TM 0 SURFACE WAVE DIFFRACTION BY A REACTIVE SURFACE WHICH 
IS TERMINATED IN A PERFECTLY-CONDUCTING SURFACE 


This chapter deals with the solution to the canonical problem 
illustrated in Fig. 3. As noted in the previous chapter, the ' • 
boundary value problems described by Figs. 2 and 4 are approxi- 
mated by the boundary value problem indicated in Fig. 3 where 
the grounded dielectric cover is replaced by a reactive surface; 
one expects such an approximation to be valid when the dielectric 
cover is sufficiently thin. This technique of replacing the . 

original boundary conditions by the approximate impedance boundary ; 
conditions for the purpose of simplifying the analysis is not new; 
it has been employed extensively in the past for analyzing cor- 
rugated surfaces, surfaces with finite conductivity, and for thin, 
grounded, dielectric slabs (a discussion on the impedance surface 
approximation for the grounded dielectric slab is available in 
Reference [8]). The main advantage of utilizing the impedance 
surface is that it simplifies the original problem to an approximate 
problem which can be solved exactly via function theoretic methods;;: 
the exact solution for the original problem is much more difficult 
as will be seen. The impedance surface for this problem is in- 
ductively reactive since the thin dielectric cover is assumed to be 
lossless. The normal surface impedance, z s associated with the 
impedance boundary (for x=0, z<0) for an assumed e-^t time dependence 
is 

(,) z s ■ - 1 x s 

The surface reactance, x g is given by 

(2 > x s = z o k x (tan \ b) Uv k d>'' 

where Zq is the free space impedance and e r is the relative 
permittivity of the dielectric slab. .The quantity k x is found 
from the TM surface wave dispersion relation given by 

k x tan k x b = e r ^(e r +l )k 2 -k x . 

The quantities k and b are the free space wave number and the thick- 
ness of the dielectric slab, respectively. Finally, 

k d =77 r k in (2). When k x b « 1, x g % z Q ( ( e r -l )/e r )kb. 

Only a single component of the magnetic field intensity exists 
for the TM excitation; for our problem, this component of the 
magnetic field is y-directed, and will be denoted by u(x,z). Let 
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u 1 (x,z) denote the y component of the magnetic field intensity of 
the incident surface wave; this surface wave is incident from the 
left as indicated in Fig. 3. In the analysis, u 1 (x,z) is allowed 
to exist even for z>0 as the "unperturbed" incident field so 
that the y component of the scattered magnetic field intensity, u s (x,z> 
accounts for the effects of the perturbation arising from the surface 
impedance discontinuity. at z=0, and x=0. The total field is thus given 
by a superposition of uMxjZ) and u s (x,z), or 

(3) u(x,z) = u^x.z) + u s (x,z). 

The form Of the incident surface wave field u 1 (x,z) is known: 

(4) u n (x,z) = e“ aX e 1 62 , for |z|<® and 0<x«» 

where 

(5) «=i« 0 z s = K 0 x S ’ 
and 

(6) e = (a 2 + k 2 ) 1/2 

with 

(7) k = to /y.e = free space wave number (y 0 and eg are the free 

J space permeability and permittivity, respec- 

tively, and uj is the angular frequency 2irf). 

a and e are the usual surface wave attenuation and propagation con- 
stants, respectively. The total field, u(x,z) satisfies the two di- 
mensional (2-D) reduced, scaler wave equation 
2 2 

(8) (— —k + + k 2 )u(x,z) = 0, 

3X 8Z 

and the following boundary conditions 

(9) 1^- + au = 0; for x=0, and z<0, 

oX 

and 

(10) It = °» for x=0 > ancl z> 0. 

oX 

In addition, (4) gives 

(11) fx _+a u i = °> for x=0 and ]z| < °° 
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and one requires that 


c ■ 0 0 1 /9 

(12) u s satisfy the radiation condition as |(x^+z z ) ' |-*» for ari^ 

e -ia)t dependent # .• 

One next defines the following Fourier transformations 

I 

isz 


(13) u(x,s) = * 


J2 7 i 


u(x,z)e dz = u + (x,s) + u_(x,s) 


where 
(14a, 14b) 


u.(x,s) e p u(x,z)e 1sz dz; u (x,s) e -~ 
jr* Jo " yf 2n 


u(x,z)e lsz dz. 


(15) u s (x,s) e -~L=- [ u s (x,z)e lsz dz = u+(x,s) + u^(x,s). 


js ; J- 


with 


(16a, 16b) 
u+(x,s) e 


u b (x,z)e lsz dz; u :> (x,s) e 


rO 


fa J- 


u s (x,z)e lsz dz. 


and 


(V7) o:(x ’ 5) s jk> 


1 r ..i/.. _\_isz _ i e' aX 


u (x,z)e dz = 

0 v 2 tt(s+s) 


(one uses (4) in obtaining the R.H.S, of (17)). The incident field, 
u 1 ’ satisfies the 2-D reduced wave equation so it follows from (8) 
that 

(18) + + k^)u s (x,z) = 0. 

ax az 

Fourier transforming (18) and using (15) one obtains 
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09 ) 


— jr + A y )u s (x,s) - 0 
dX 

where 

(20) J7 =Jk 2 -s 2 . 

A > 


The general solution to (19) is given by 


(21) 


u s (x,s) = 


aC aC 

u + + u _ 


iJT x 
= C(s)e x 


-i Jk x 

D(s)e x 


In order to simplify some of the analysis, one introduces a small 
loss in the medium above x=0; i.e., one lets k (and x x ) be complex 

(22) k = k-j + i k 2> (k 2 << k^ ; k^ > 0, k 2 > 0). 

This loss may be set equal to zero once the analysis is completed. 

From (12) one notes that u s must satisfy the radiation condition 
as x oo so that D(s) = 0. Thus 

i/T v x 

(23) u s (x,s) = C(s)e x 
wi th 

(24) =J k 2 -s 2 = +i J s 2 -k 2 , 

where the choice. of the branch of JY X is the one for which Im JY^ > 0. 
One notes that u 1 (x,z) already satisfies the impedance boundary con- 
dition; therefore 

(25) + a u s = 0; for x = 0, and z < 0 . 

o X 

Furthermore, (10) indicates that 

(26) ~~ = - ; for x = 0, and z > 0 . 
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In order to apply (25) and (26) to u s (x,z), one takes the inverse ^ 
Fourier transform of (23) to obtain 


u s (x,z)=-M u s (x,s)e~ 1sz 

-• ‘ • v/2* ” 


ds 


or 


(27) u s (x,z) = 


J 2tt 


i/Tx . ; . ;; * 

C(s)e x e" lsz ds 1 ' 


i' S' \ 


> r ir,\ 


:f.Fc 


•ill •/ 


Substituting the result on the R.H.S. of (27) into. (25) V; . 
and (26) leads to the following pair of integral equations 


(28) 


Q\ -i«) C(s)e" isz ds = 0, z<0 

A . ,* 


(29) f ^ c (s)e“ isz ds = -ic«/2^ e iBZ , z>0 

It is convenient to define 

(30) <t>(s) =JT X C(s) 

and 







so that an alternate set of dual integral equations is obtained by 
substituting (30) and (31) into (28) and (29); thus 


(32) 


*00 1 • 

L" (s) <|>(s)e“ lsz ds = 0, z < 0 

— 00 



(33) 


<j>(s)e" 1sz ds = -i Jzn ae i0z , z > 0 


12 



One may now solve for <f>(s ) given in (32) and (33) via a factorization 
procedure. 


Let 


(34) s = o+iT, 

The region in the complex s plane for t > -k 2 will be denoted by 
U.H.P. (upper half s plane), and the region in the complex s 
plane for t < k2 will be denoted by L.H.P. (lower half s plane). 
Initially, one assumes that L _1 (s) 4 > (s ) -* 0 in the U.H.P. as 
| s | -*■ ® there, and one also requires that <j» (s ) -*• 0 in the. L.H.P. 
as | s | -*■ » there. It will be indicated aposteriori that these 
assumptions are true. Under these assumptions, one may employ 
Jordon's lemma to close the contour of integration in, the U.H.P. 
for (32) and similarly close the contour of integration in the 
L.H.P. for (33) so that 

L -1 (s) <j>(s)e" isz ds = 0, z < 0 
J 

e> 

(36) j <f>(s)e -lsz ds = -i|2ir ae 10z , z > 0 

$ 

From (36), one requires that <j>(s) have a pole in the L.H.P. to 
provide a residue contribution corresponding to cte 1 "# 2 . From 

(35) one may write 

(37) l-’(s) *(s) = A + (s) 

where A*(s) is as yet an unknown analytic function in the U.H.P. 
Similarly, one may define another as yet unknown analytic function, 
A_(s) in the L.H.P. such that 

A_(s) 

^ JTH A_(-B)(s+e) = 

If one employs (38) in the L.H.S. of (36), one then obtains the 
value on the R.H.S. of (36) via the Cauchy Residue Theorem. 
Eliminating <j>(s) between (37) and (38) yields 


(35) 

and 
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m L(s) "jf; TOT W w 

The unknowns Af(s) and A_(s) can be found by factorizing L(s) into 

L+(s) I (s ) (where L+(s) is analytic in the U.H.P. and 1 (s ) is 

analytic in the L.H.P.; L+(s) are referred to as the Wiener-Hopf 
factors). Clearly I (s ) is proportional to A_(s)/A_(-b); hence 

L_(s) A (s) 

(40) 1TF5T' TTFbT , 

(The proportionality constant in the above equation is unity since 
the R.H.S. equals unity where s = -B ). Thus, it follows that 


(41) A (s) = — 2 

jTjrL_(-6)L + (s) (s+b) 

The explicit forms of the functions L+(s) and I ( s ) are derived in 

Section II of Appendix I. These L ± (s) functions are deduced from 
an appropriate limiting procedure on the factorization function 
for a related closed region boundary value problem; this technique 
leads to factors which are convenient for numerical processing. One 
notes that if the boundary conditions (25) and (26) were employed in 
the transformed domain (x,s-domain) instead of in the spatial domain 
(x,z-domain) as done in this analysis, one would have obtained a . 
Wiener-Hopf equation for the unknown <j>(s). However, the solution 
to the Wiener-Hopf equation requires decomposition functions in 
addition to factorization; thus, by formulating this problem in 
terms of the dual integral equations one avoids any decomposition. 
The solution to the Wiener-Hopf equation also requires an analytic 
continuation argument which is not essential for the solution of 
dual integral equations by the factorization procedure. Some 
author's prefer to use the analytic continuation argument even to 
solve dual integral equations by the factorization method; in 
contrast, the procedure outlined above allows one to by-pass the 
analytic continuation argument to solve for A + (s). Incorporating 

(41) in (30) via (37) yields 

L (s) 

(42) C(s) =-$=- _L 

^‘■.(-e) 777 2 (s +8 ) • 

One may next utilize (42) in (23) to obtain 


14 



( 43 ) 


^ e^i 


2 2 
k -s x 


^ L >»> J?7 2 ( s +6 ) 


Inverse Fourier transforming (43) leads to a formal solution for 
u s (x,z) in, terms of an inversion integral with known terms in 
the integrand; 


(44) u (x,z) = 27 


, r L.(s)e 1 J k2 - 


s 2 x 


L_(-6)Jk 2 -s 2 (s+B) 


-isz . 

e ds . . 


An evaluation of the integral in (44) leads to the complete solution 
of our problem. Before procedi ng further, one may check aposteriori 
if the assumptions concerning the behavior of 4>(s) in the L.H.P., 
and of L“ ' (s)<j>(s) in the U.H.P., have all been satisfied. From the 
asymptotic behavior of the factors L ± (|s|-*=°) as given in Appendix I 
(Section II) - one notes that L+(s) -v (constant). Hence, from (41), 
A+(s) 'O 0(s ■“'■)* so that L" 1 (s)(f>(s) ^0(s -1 ) and <j>(s) 'x-0(s - '). It 
is now easily verified from these asymptotic behaviors that all the 
previous assumptions are indeed true. The contour of integration 
in (44) may be closed in the U.H.P., and the Cauchy residue theorem 
applied to obtain the fields which exist in the region z < 0, x _> 0; - 
whereas one may close the contour, of integration in the L.H.P. to 
obtain fields which exist in the region z > 0, and x 0. The 
singularities of the integrand of (44) in the complex s plane and 
the deformed paths of integration are indicated in Fig. 5. One 

notes that 1 (s ) contributes the pole at s = -3 (see (44)). The 

original contour may be deformed to one in the L.H.P. as follows 


(45) 



/•CO 




= - ;■ 


+ 

V -J 

— 00 

. 

C z>0 1 


C br. 1 +C br. 2 



for z > 0 


Similarly, 


the deformation into the U.H.P. is given by 


(46) 


•CO 


+ f 

J — oo 

J 

o 

N 

A 

o 


'br. 1 


+ C' 
u br. 2 



for z < 0. 


The contributions to the integral from 
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Fig. 5 —Contours of integration in the complex s-plane 
which occur in the solution to the problem in 

Fig. 3. 



vanish. The integrals over Ct> r . 1 + Cfc>r. 2 and c br. l+C^r. ? yield 
the radiation field which may be asymptotically evaluated in closed, 
form yi a the saddle point method; whereas, the integrals over Cg 
and Cg yield the transmitted and reflected surface waves, respectively. 
There can be no transmi tted surface wave above the perfectly con- 
ducting surface (x=0, z>0); however, the transmitted surface wave 
arising from the integration over the circular contour C„ exactly 
cancels the incident surface wave field which was assumed to exist 
for z>0 in the analysis. The transmitted surface wave field for 
z>0 is given by the residue at s— $ as: . • . 
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( 47 ) 


2n 



aL (s)e 


W 


k 2 -s 2 x 


L_ {— B W k 2 -s 2 (s+e ) 



-2771 


«L_(-e) 

L ( - g ) i a 2n e 


- aX e 1Bz =-u 


i 


J 2 2 

k -g = ia has been used in deriving the above result) 
Similarly, the reflected surface wave for z<0 is given by the residue 
at s=g as: 


(48) 


h 

L 


L (s)e 




s 2 x 


J-e)J 


k 2 -s 2 (s+g) 


-isz . <(6)e“ ax _ i62 

6 dS 2lT V(-B)ia(2Tr)(23) 6 


where 



from Appendix I (Section II). The symbol P.V. in front of the integral 
denotes the principal value of the integral. The quantity £ (g) may 
be identified as 


(50) £(g) = - Residue of L_(g) = - lim (s-g) L (s) 

s+g 



Since the value of the surface wave field at x=0, z=0 is unity, the 
reflected surface wave denoted by u r (x,z) for z <0 may be expressed 
as 


(52) 


u r (x,z) = R 1 


-aX 



» 
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where Ri denotes the surface wave reflection coefficient for this 
canonical problem. From (48), (49) and (50), R-, is explicitly 
given by 


(53) 






2 2 

Nk -q +e- 


J~2 

q +a 


dq . 


The radiation fields obtainable from C^ r ] + Ck r> i and c h r l + 

Cbr. 2 via an asymptotic evaluation are the fields diffracted by the 
impedance discontinuity at (x=0, z=0). It is sufficient -.to evaluate 
the branch cut integral either in the U.H.P. or in the L.H.P. by the 
saddle point method (both lead to the same result). Instead of 
asymptotically evaluating the integrals by deforming the contours 
Cbr. 1 + ^br. 2 into a steepest descent contour which passes through 
the ’saddle point in the s-plane, it appears to be more convenient to 
employ a polar transformation in the integral of (44) and later evaluate 
the transformed integral via the saddle point method. The use of polar 
transformations is a common procedure for evaluating integrals of the 
type given in (44) (see [8,9] for example). The polar transformations 
are V" 


(54) s = k sin 5;v s 2 -k 2 = -ik cos 5 


(55) x = p cos 4> and z = p sin <j> . 

Figure 6 illustrates the quantities p and <j>; these correspond to the 
polar coordinates at the field point. The polar transformation yields 
a single valued integrand in the transformed plane (complex 5 -plane). 
The new contour in the 5 -plane may be deformed into the steepest 
descent path denoted by Csgp in order to obtain the saddle point 
contribution to the integral. The saddle point contribution gives the 
diffracted field contribution; whereas, the residues corresponding to 
the poles; which may be crossed in deforming the original contour in the 
5 -plane to the G$dp contour give rise to surface wave contributions to 
the field. In this problem, these pole contributions in the 5 -plane 
are identical to the pole contributions arising from contours Cg 
and Cg in the s-plane as given by (47) and (48); hence, these will 
not be calculated again. The complex s-plane representation appears 
more suitable for calculating surface wave contributions to the 
field; whereas, the 5 -plane representation appears to be more 
convenient for calculating the diffracted field. Thedetails of 
the transformation of the integral in (44) from the s-plane to the 
5 -plane, and the complex 5 -plane topology will not be presented 
for the sake of brevity. The final expression for the integral over 
the steepest descent path is denoted by ud, where 


18 




Fig. 6 — Polar coordinates for the field point. 


(56) 


u d (p,*) 



L_ (k sin ?) 

L C-eMe+k sin 5 ) 


e ikp cos( 5 +*) d 


The above integral is evaluated asymptotically by the saddle point 
method; and the leading term in the resulting asymptotic solution 
for large kp is 


(57) 


U d (p,*) 



JJ^k 


L_(-k sin tj>) e ikp 
L_(- 3 ) (B-k sin <j>) * 


Since the value of the incident surface wave field at x=0, and z=0 
is unity, one may express u^ as 

j -ikp _ . ; 

(58) u d p.D 1 (»)^ 

Jp 


where D-j represents the surface wave diffraction coefficient for 
this canonical problem. From (57) and (58) D-j is thus given by 


(59) 


-i - 
1 4 


ae • L_(-k S1 " n 4 >) 

D l^ = L .(-B)(B-k sin 


Substituting for L_(s), D-j ( <}>) becomes 
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CHAPTER III 

TM 0 SURFACE WAVE DIFFRACTION BY A RECESSED SEMI-INFINITE 
DIELECTRIC SLAB WHICH IS FLUSH MOUNTED IN 
A PERFECTLY CONDUCTING SURFACE 


This chapter deals with the solution to the canonical problem 
illustrated in Fig. 2. Two approaches for solving this boundary 
value problem are explored. One approach is based on the generalized 
scattering matrix technique (henceforth denoted by GSMT) which was., 
recently developed by Pace and Mittra[7] for analysing a class of 
waveguide discontinuity problems; the generalized scattering 
matrix approach leads to a formally exact solution for the problem. 
The scattering matrices involved are found from the solution of 
associated problems by the Wiener-Hopf method. The second approach 
is based on the variational technique which leads to an approximate 
solution. The details pertaining to the solution of this canonical 
problem via the GSMT is indicated in section A of this chapter; 
whereas, section B of this chapter describes the variational 
solution to the same problem. 

A. Solution Based on the GSMT 

The GSMT was originally developed for analyzing the effects of 
discontinuities in closed waveguide regions [7]; however, more 
recently, Lee and Mittra[10] extended the application of this 
technique by solving an open region problem involving the dif- 
fraction of a plane wave by a thick conducting half plane (they 
actually considered plane wave diffraction by a dielectric loaded 
parallel plate waveguide which reduces to the case of a thick 
conducting half plane when the relative dielectric constant tends 
to infinity). Kashyap and Hamid [11] later extended the results 
of Lee and Mittra to treat the diffraction by a slit in a thick- 
conducting screen via the GSMT. The concept of the generalized 
scattering matrix is closely connected with the scattering matrix 
theory used for describing microwave networks. The basic dif- 
ference between the ordinary scattering matrix and the generalized 
scattering matrix is that the former considers only propagating 
modes; whereas, the latter considers both evanescent as well as 
propagating modes. As a result of the extension to include 
evanescent modes, the generalized scattering matrix is of infinite 
order, and is non-symmetric (the symmetry in the ordinary 
scattering matrix results from the normalization on the propa- 
gating modes such that they carry unit power; such a normalization 
cannot be applied to evanescent modes). In general, boundary value 
problems whose geometrical configurations can be identified as 
being composed of two or more junctions can be solved via the GSMT, 
provided that the generalized scattering matrices associated with 
these junctions can be found. In the GSMT, the formally exact 
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solution to the boundary value problem is obtained by systematically 
summing up field contributions arising from the multiple scattering 
between these junctions. The multiple scattering process between the-, 
junctions can be described by a convergent Neumann series whose terms 
contain these generalized scattering matrices [7]; however, one can 
directly arrive at the same result via a self consistent method for , 
multiple interactions. The latter viewpoint will be employed in our v . 
analysis. The formal solution arising from the multiple interactions;; 
contains infinite order matrices. In most practical cases it is 
possible to obtain accurate results by truncating the matrices ! 

to a reasonably small size as seen from [10] and from the results 
obtained in this report. V.' 

At first glance, the canonical problem illustrated in Fig. 2 does 
not appear to be composed of two junctions. However, if one assumes 
that the dielectric slab is extended by an "infinitesimal" distance,.. 

6 into the conducting region (i.e., into the region 0 <_z < 6; 

0 < x <^b), one creates physically identifiable junctions "© and 
(2J* as indicated in Fig. 7. When 6 0, the problem illustrated 

in Fig. 7 becomes identical to the canonical problem illustrated in 
Fig. 2. For later convenience, one defines the regions ®, (|) 
and © surrounding junction (l) such that region ® exists for 
(x > b; Izj < «), region ® exists for (0 < x < b; z < 0), and 
region (t) exists for (0 < x < b; 0 < z < <S). These regions are also 
indicated in Fig. 7. The geometries of junctions ©and (2) are in- 
dicated indicated in Figs. 8a and 8b, respectively. The properties 
of junctions 0 and © may now be described in terms of their:, 
respective generalized scattering matrices. 



Fig. 7 —Junctions © and © for the recessed slab 
in Fig. 2. 
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In order to introduce the generalized scattering matrices 
for junction 0, one begins by considering the geometry for this - 
junction alone as shown in Fig. 8a. Junction 0 is, excited by the 
TM 0 surface wave which is incident from the left as in Fig. 7. This 
incident surface wave is scattered by 0 such that a part of it is 
scattered into region 0 and the rest is scattered into regions . ®. 
and ®. The fields transmitted into region (C) are described by. the 
generalized scattering matrix denoted by S2i; whereas, the fields 
scattered into regions ® and (§) may be described by the • 
generalized scattering matrix denoted by Sy|. The scattering , 
processes described by S2i and Sn are indicated in Fig. 9a. The . 
field transmitted into region © impinges on junction © from which 
it is reflected towards ® . The field in region © is represented 
by a sum of TM on waveguide modes (for n = 0, 1, 2, ••• °°). Thus,, 
any TMon mode which is incident at 0 after reflection at © is 
partly scattered back into the waveguide region © , and the rest is 
scattered into regions ® and ® . The field scattered back into © 
is described by the generalized scattering matrix denoted by S2?; 
whereas, the corresponding scattered fields in regions ® and © are 
described by the generalized scattering matrix denoted by S-j 2- The 
scattering processes described by S22 and Si 2 are shown in Fig. 9b. 

It is clear that the matrices Sn and S21 may be obtained from the 
solution to the problem illustrated in Fig. 9a; whereas, the matrices 
S22 and S12 may be obtained from the solution to the problem illustrated 
in Fig. 9b. , The problems illustrated in Figs. 9a and 9b may be re- 
garded as auxiliary problems associated with the original canonical- 
problem; the solutions to these auxiliary canonical ; problems are 
instrumental to the solution of the original canonical problem 
illustrated in Fig. 2. The generalized scattering matrix for junction 
® is denoted by r.< r corresponds to the reflection coefficients of 
the TM on modes which are reflected at ©; it is easily seen that r 
is an identity matrix of infinite order. 

The explicit forms of the generalized scattering matrices 
Sn, S12* S21 and S22 are derived next. The derivations for Sn 
and S21 are indicated first; these are followed by the derivations 
for Si 2 and S22. 'It is a trivial matter to determine r, and hence 
the derivation for r is not presented here for the sake of brevity. 

The surface wave reflection' and diffraction coefficients are then 
constructed from these scattering matrices by considering inter- 
actions between the junctions 0 and ©. 

Derivation of S^ and S 2 -| : 

The generalized scattering matrices Sn and S21 are obtained 
from the solution to the auxiliary problem illustrated in 
Fig. 9a as mentioned before. The excitation for this problem is a 
TM q surface wave incident at 0 from the left as^shown in Fig. 7. 

For the TM case, the magnetic field intensity is y-directed; 


23 




(a) (b! 

Fig. 8 — Geometries associated with the individual junctions. 



(INTO REGION ©) 


INCIDENT FIELDH 


,-a (x-b) 4- !/?z 


,x>b 


COS k** ; a. 


S„> 

(INTO REGION @ 
AND ® ) 


(a) 



INCIDENT FIELD = COS ~ xe -l ^ nZ 

D 




(INTO REGIONS (INTO REGION ©) 
(A) AND ® ) 


(b) 


Fig. 9 — Auxiliary canonical problems associated with 
the recessed slab problem. 
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hence, in the following analysis one defines the scaler fields , 
<j>, and to respectively correspond to the y-directed components 
of the incident, scattered and the total magnetic field intensities 
For the purposes of anlaysis, it is convenient to assume that 
exists even for z > 0. The final solution will contain a con- 
tribution to o which exactly cancels the incident field for z > 0. 
The form of the incident field s 1 is 


(61) 


with 

(62) 


f 


e -a(x-b)+i gz 



cos k x x J gz 
cos k x b e 


x >_ b, |z | < oo 

0<_X£b, |z| < co 





(63) a 2 + k 2 = g 2 , 


where a and g are the TM 0 surface wave attenuation and propagation 
constants which are found from the TM surface wave dispersion 
relation 


(64) k x tan k y b = e r a 

k and kd are the wave numbers in the regions x > b, and 0 < x < b, 
respectively. Thus, 


(65) k^ J!~ k, 

where e r is the relative permittivity of the dielectric in 0<x<b. 
The total field is a superposition of $>i and $. 

( 66 ) ^ + $. 

o 1 , and $ satisfy the homogeneous, 2-D, reduced wave equation. 
2 2 \ 

(67) + ^-y+ k 2 ]$(x,z) = 0; x > b, |z| < - 

V ax az J 
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and 


■ w 


(68) 


( 


2 2 \ 

+ -^-k- . + J $(x,z) = 0; 0<_x 7 £b, 

ax az ■/ ' ;; i • 


In addition, $ must satisfy the following; boundary conditions and 
the radiation condition. 


< 69 > If 


x=0 


=0; for | z | <oo 


(70) fj {»+ t 1 } | x=b =0 = 1- |j{» + f 1 } l x=b for z>0 

+ r 


( 71 ) $ lx=b + =$ lx=b _ ;for z< ? 


(72) 


_a$ 

ax 


_ 1 _ 21 
x=b + " e r 9X 


x=b 


; for z < 0 


2 2 1/2 

(73) $ satisfies. the radiation condition as |(x +z ) | -> 

for an e - "*^ time dependence. 

One defines the following Fourier transformations:,.! '■ 


1 


(74) i(x,s) = J $(x,z)e lsz dz - $ + (x,s) +: i_(x,s) 


where 
(75a, 75b) 


.to # -yj # 

3 (x,s) 5 • $(x,z)e 1sz dz; 3_(x,s) e — Lr $(x,z)e lsz dz 

J 2tT 0 . >/2 7T , 00 
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and s = a + ir as in Chapter II. Also 
(76a, 76b) 

$+(x,s) = -=f $ 1 (x,z)e 1sz dz; V (x,s) = 

j2Jo 


JTn i 


r oo 1 

3X 


a -*- e 1sz dz 


From (61) and (76a), one easily obtains 
r 


and 


(77) $’ = 


e -a(x~ b ) 

s/i7 ' 5+6 

1 cos k x* i 

nr cos kb s +3 * 

yJCVi X 


r 


a e -a(x-b) _i_ 


v , & ~ 5+6 * 

< 78 > *1 H k sin k x . 




X_ 1 

cos kb s+6* 


x > b 


0 < x < b 


x > b 


0 < x < b 


Fourier transforming (67) and (68) w.r.t. ? yields 


(79) 
and 

(80) 


i 

1 


4 - (s 2 -k 2 ) 
ax 


f 3 , (X > ,S) = 0 


.2. 


ax 


- (= 2 -^ d ) 


$(x < ,s) = 0 


Where x> and x< denote values of x for x >_ b and 0 x £ b, 
respectively; thus, x> ->• b + and x < b. as x b. General 
solutions to (79) and (80) are 


(81) S( X> ,s) = A(s)e" YX + B(s)e yx 

(82) i(x < ,s) = C(s) cosh y-|X + D(s) sinh y-|X 
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where 


(83) y = js 2 - k 2 = - i Jk*~- s 2 . . . , 

and 

(84) y, 5 Js 2 - k 2 

' ;> \ 

Here, the choice of the branch of y is the one for which Re y > 0.' 
For ease of analysis, one introduces a small loss into the regions 
x ^ b and 0 £ x ^ b; i.e., one lets k and be complex. 

(85) k = k-j + ik 2 , (k 2 << k-j ; k^ > 0, k 2 > 0) 
and 

(86) k d = k 3 + ik 4 , (k 4 « k 3 ; k 3 >0, k 4 > 0) 

The radiation condition (73) requires that B(s) = 0; and (69) 
requires that D(s) = 0. Thus, . v . 

(87) $( X> ) = A(s)e" YX = * + ( Xj j + S_( X> ) 

(88) *(x < ) = C(s) cosh y-j x = S + (x < ) + $_(x < ) 

Employing (71) in (87) and (88) yields 

(89) 3 + (b + ) " = A (s)e“ Yb - C(s) cosh y^b 

Similarly, from (72), (87) and (88) one obtains 

(90) e r $J.( b + ) -*|(b_) = -e r YA(s)e" Yb -Y-|C(s)sinh y-jb . 
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where primes indicate differentiation w.r.t. x. Employing (70), 
and (78) in (90), one arrives at the following result after some 
manipulations (and the use of (64)) 

(91) £r yA(s)e‘ Yb = ~Y-j C(s) sinh Y-jb 


One next defines 


(92) J + (s) = 5 + (bJ - 3 + (b + ) 

Utilizing (91) and (92) into (89) gives 

(93) -J + (s) - YA(s)e _Tb G _1 (s ) 
where 

(94) G(s) = (yy-, sinh Y-jb)^ sinh Y^b + e r Y cosh Y-jb)” 1 
From the fact that 

(95) $;(b + ) + *:(b + ) = -YA(s)e" Yb , 
and 

(96) ;;<b + ) « -*j'(b + ) , 
one obtains 

(97) -«{' (b + ) + $^(b + ) » -vA(s)c' yb . 

Finally, incorporating (97) into (93) and using (78) allows one to 
obtain a Wiener-Hopf equation 
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(98) J. (s') G(s ) =. »>. ) + -j? 

fa (s+e) ' . , 

From the behavior of <|>(x,z) and 3.<()(x,z)/3x as |z [ -* ”, it fs evident 
that the transform $(x>) is analytic in |x| < k 2 *, whereas, $(x<) is 
analytic in |x| < (< 4 . Also, J+(s) is analytic for x > -[Min(k 2 ,k 4 )]. 

G (s ) is analytic for |x| <Min(k 2 ,k 4 ), and $i(b+) is analytic for x<k 2 . 
Thus, (98) is analytic in the strip ft | < Min(k 2 ,k 4 ). One begins 
solving (98) by factorizing G(s), i.e. -,i 


(99) G(s) = G + (s) G_(s) 


Where G+ is analytic in the U.H.P. (i.e., for x > - Min(k 2 ,k 4 )) and 
G_ is analytic in the L.H.P. (i.e, for x < Min(k 2 ,k 4 )). The 
factorization functions G±(s) for this problem are derived in 
Appendix II via a formal factorization procedure since the limiting 
type factorization scheme indicated in Appendix I is shown' to cause 
convergence problems in the integral representations for the factors 
G+(s) (see section III of Appendix I for these details). The use 
of (99) in (98) yields 


( 100 ) 


*:o>+> 

G + (s)J + (s) = £j s -) + 


ict 

^G_(s)(s+3) 


The R.H.S. of (100) is analytic in L.H.P. except at s=- 3 . It is 
convenient to put the preceeding equation into the form 


001) 

G + (s)J + (s) 


_ ig 

J27g_(-6)(s+B) 


G - (S) _J27(s+6) 


(g; 1 (s)-g ; 1 (-e )3 


The L.H.S. of (101) is analytic in the U.H.P. ; whereas, the R.H.S. 
of (101) is now analytic in the L.H.P. Hence, the R.H.S. and the ; 
L.H.S. of (101) together define an entire function denoted by P(s). 
Since G(s) ^ 0(s) as |s|-*» in the strip of analyticity, G ± (s) % O(s^) 
as [ s | », in their regions of analyticity. The edge condition for 

this problem dictates that the magnetic field intensity must behave 
as p + ^/2 -j n neighborhood of the edge as p -*• 0 (p = radial distance 
away from the edge at x = 0, z = 0), which correspondingly implies 
that J+(s) ^ 0(s“3/2) as ] s | « in the strip (since 0+(s) may be 

interpreted as the transform of the current density on the plate at 
x = b for z > 0); the edge condition for this problem is deduced by 
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an application of the technique described in [9], The entire function 
P(s) therefore equals zero from the asymptotic conditions for large s. 
Hence, each side of (101) is zero. It follows that 

(102) J + (s) = — -12 

>^n"G + (s)G_(-8) (s+e) 

Utilizing (102) in (93) gives 
-iae Yb G (s) 

(103) A(s) = — 

>/2^YG_(-g)(s+g) 


From (103) and (91 ), one can obtain an expression for C(s) in 
terms of known quantities. Incorporating this value of C(s) in (87), 
and the above value of A(s) in (88) one obtains 

-iae Yb G (s)e' YX > 

(104) $(x.s) = — , 

J2^y G ( — 0 ) ( s +$) 


and 


( 1 05 )r $ ( x < , s ) 


1e r a G _( s )cosh y 1 x < 
^rY-|SinhY^b G_(-g)(s+g) 


Finally, inverse transforming (104) and (105) yields 


, -i a G (s)e 
(106) $( X> ,z) = ~ 1 


~y(x. - b) 


-a. yG (-g)(s+g) 


e~ isz ds 


and 


•, f - ie aG (s)cosh Y,X^ 
(107) *(x <t z) f r - 1 < 


-co Y-|Sinh y-j b G_(-e)(s+e) 


„-isz 
e ds 
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thereby completing the formal solution to this canonical problem. 

One may close the integration contours along the real axis in (106) 
and (107) by semi-circular contours at infinity (except for the 
deformations around the branch cuts) so that the fields represented 
by <i>(x > ,z) and $(x < ,z) may be calculated via Cauchy's residue 
theorem. The semi-circular contours at infinity provide vanishing 
contributions to the integral. If one closes these contours 
defined along the real axis by a semi-circular contour at infinity 
in the U.H.P., the fields for the region z < 0 are obtained; whereas 
if one closes the path of integration defined on the real axis by a 
semi-circular contour at infinity in the L.H.P., the fields for 
the region z > 0 are obtained. The pertinent contours in the 
complex s-plane are shown in Fig. 10. Clearly, the residue at the 
pole s = -6 in the L.H.P. which is obtained by closing the contour 
of integration in the L.H.P. exactly cancels the incident field 
for z > 0, as was anticipated. 



Fig. 10 — Contours of integration in the complex s-plane 
which occur in the solution to the auxiliary 
canon 'cal problems. 
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The diffracted field in region (§) may be obtained from (106) 
via an asymptotic approximation to the integral. This diffracted 
field arises from the branch cut contributions to (106). One 
begins by introducing the usual polar transformations 


(108a) 


s = k sin a. 



= k COS a 


and 


(108b) x = p cos <j>. 


z = p sin 4 > . 


Incorporating the above transformations into (106) leads to the 
following integral along the contour (> in the complex g-plane 


(109) *(p,<fr) 


a 

2tt 


G (k sin a)e 


ikpcos(a+<f>) 


•’(> G_(-g)(k sin a + g) 


-i kb cos a.- 
e da 


The quantities p and <p are indicated in Fig. 8a. The contour O 
may be deformed into the steepest descent contour for the above 0 
integral in the a-plane, and the resulting integral may then be 
asymptotically evaluated for large kp by the saddle point method, 
to give the result for the diffracted field denoted by $ d as 


(no) 


4> d (p,<f>) 




f\*) 


e ikp 

7T 


where 


(HI) 


, aG (-k sin 4> )e" lkb cos * -i l 

F A U) = e 4 . 

>/2TTk G_(-g)(g-k sin <j>) 


Since the surface wave incident on the discontinuity at z = 0 in 
Fig. 8a has unit amplitude, it is apparent that FA^) is a diffraction 
coefficient for this auxiliary problem. It is interesting that it has 
the same form as the diffraction coefficient Di for the diffraction of 
a surface wave by a discontinuity in surface reactance (see (59)) 
except for the factor e-ikb cos which is due to a difference in 
the origin of the coordinates for the two problems. 
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The TM 0 surface wave which is reflected from 0 is obtained 
from the residue at the pole s = s of G_(s) in (106) for the 
reflected surface wave in the region x ^ b and z < 0 (i.e., for 
the part of region ® for which z < 0), or from the residue at 
the pole s=3 of G_(s) in (107) for the reflected surface wave 
in region® (i.e., 0 x < b; z < 0). The result for the r 

reflected surface wave field in region 0 (for z < 0) which is 
denoted by $ r is indicated below; it is obtained in a straight- 
forward manner by dosing the contour of integration in (106) by ' 1 
a semi-circular contour at infinity in the U.H.P. and evaluating 
the residue at s = 3, and by replacing G (s) in (106) by G(s)/G + (s) 
for convenience. 


( 112 ) 


Tg e -a(x-b) e' Uz 
' t = G_(-$)G + (s)^6 


x > b, z < 0 


where 


(113a) ^ = Residue of G($) = lim (s-b) G(s) 

° s->3 

Using (94), is found to be 


(113b) = 


a(ajl 


“ 0 2 ) 2 


e[e r (kj-k 2 )+ab(k^-B 2 +e^a 2 )] 


One may now express v above as 


(114) = R s « e -«(x- b ) e -4SZ 


; x > b, z < 0 


where R sw is the TM Q surface wave reflection coefficient at 0 , 


(115) R sw = 


7 ? 


% 


2bG_(-3)G + (b) 2b g;(b) 
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(since G_(-s) = G + (s) as indicated in Appendix II). 

The fields transmitted into region © are given by the 
residues at the poles of ( yj sinh Ylb)“‘ in the integrand of 
(107). Before evaluating this contribution it is convenient to 
re-write (107) as 


-1 


(116) $(x < ,z) = ^ 


00 i e„aG 


_(s)cosJk^- 


2 -isz 
s x,, e 


~ m J k^-s 2 sinjk 2 -s 2 b G + (g)(s+3) 


ds 


f ~ p ? 12 2 

The zeros of'Jkj-s sin^k^-s b corresponding to the fields for z>0 


occur at s =-3 . 3„ is given by 

mm 


(117) 


3, 


m 


-k-(¥f ■■ 


m = 0, 1 , 2, 


The fields transmitted into the waveguide region 0 are thus 
obtained by closing the contour for (116) in the L.H.P. and by 
summing up the residues contributions arising from the poles of 

(>/ k d“S 2 sinjk 2 -s 2 b) -1 at s = 

Let the transmitted field in region 0 be denoted by 

i 3_z 


(118) 


00 G (-3 )cos x e 
.t _ r - m b < 

<J> - -e„a I 


m 


m=0 


{^s (l^d- 2 sin / k2 d- 2b )K G +(3)(3-B m ) 

N p m 


or 


ik d z 


(119) 9 Z =< 


e r aG +^ ( j) e 

2bk d G + (B)(f5-k d ) ’ 


for m = 0 


« G (-e )cos 52- x e 1 V 

- £ a l ^ , for rfO 

m=l (-1) b3_G + (s)(B-6 ) and m=l ,2,3,-** 


f may be expressed as 


. 1 3 z 

(120) = l T m cos j— x < e m ; for 0 <_ x <_ b and z > 0 


m=0 
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where T m denote the modal coefficients of the fields transmitted 
into region 0 . From (119) and (120) one may define T m ~as 


( 121 ) 


e aG,(g ) 
r + nr 


m 


4 n,(-')\ b G + < 6 ><S-eJ 


wi th 



for m = 0 
for m t 0 



The generalized scattering matrix sn for junction 0 
describes the fields scattered in regions (S) and 0 when a 
TM 0 surface wave is incident at junction 0; this scattered field ' 
is composed of a reflected TM 0 surface wave field and a diffracted 
field contribution. If one is interested in the reflected surface 
wave calculation, one defines s-n to be R sw of (115); on the other 
hand, if one is interested in calculating the diffracted field in 
region @ , one defines s^ to be F A (<j>) of (111). Thus, 

R sw , for reflected surface wave field calculations 
in regions ® (for z<0), and ® 

N . 

F A (<j>), for diffracted field calculations in 
region (§) . , 


(122) s-i -i = j 


The generalized scattering matrix S2l represents the modal.. _ 
coefficients of the field transmitted inio the waveguide region 0 
when a TM 0 surface wave is incident at 0 from the left. The field 
transmitted into © is given by (120). Hence, s 21 may be 
conveniently represented by a matrix whose elements are T 0 , Ti , T 2 , 
... Tm •••» where T m is defined in (121). Thus, one defines 

S21 3 S 


(123) S 21 - [T Q T^ Tg 


m 


] 


T 

IxM 


where T is the transpose operator. s ?1 is an M x 1 matrix in 
which M •* ». 
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Derivation of and s^ 


The generalized scattering matrices s-|2 and S 22 are obtained 
from the solution of the auxiliary problem illustrated in Fig. 9b. 
The excitation for this problem is a TM on waveguide mode which is 
incident at G) from the waveguide region © (i.e., from the 
right). One denotes the y-component of the incident, scattered 
and the total magnetic field intensities by $i, $ n and . 
respectively. Then <i> n = ^ where the incident field, is 

assumed to exist even for z < 0 for convenience of analysis. The 
final solution will contain a contribution to $ n which exactly 
cancels. the incident field for z < 0. The form of the incident 
field *1 is 


-l $ z 

(124) <J>^' = cos x e n ; for 0 <_ x b and |z| < <» 


where 


< 125 > e n “j k d-({r|- "- 0 * 1 - 2 "- 


and kd has been defined in (65). Bates and Mittra [2], and Angulo 
and Chang [15], have solved the problem illustrated in Fig. 9b for 
the special case when the structure is excited by a TEM waveguide 
mode; this corresponds to the n = 0 case in our problem. <*>i 
satisfies the 2-D, reduced wave equation 


(126) 



0 <_ x . <_ b, |z| <« 


= 0 everywhere else. Also, satisfies the 2-D, reduced wave 
equation; hence, $ n must satisfy 


027) 


+ - — + k 2 ) $ (x,z) = 0; for x >_ b, |z| < <*> 

V 9x 3z J n 


and 

(128) 




$ n (x,z) = 0; 


for 0<^x<^b, J z | < oo 
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In addition, ^ satisfies the following boundary conditions and 
the radiation condition. 


3$ 

029) 


x=b 


= 0, for z > 0 


(130) 


3$ 

r 

3 X 


x=0 


= 0, for |z| < « 


<131) * nl x=b, = ( ®n + ’ forz<0 


3$ n 

< 132 ) 1ST 


. = ( + $ ) , , for z < 0 

x=b. p 9 x ' n *n ;| x=b * 

+ r 


2 2 1/2 

(133) $ satisfies the radiation condition as |(x +z ) |-*» 

for an e~ 1a)t time dependence. 


As before, one defines 


(134) 3„(x,s) = 7 = I $ n (x,z)e isz dz = 8 *(x,s) + Ox,s) 


~JU - - " 


where 

(135) 

A+ / X 1 

, n ( x,s) ^ 


r co 


$ (x,z)e 1 sz dz; T(x,s) _ — . 

o n n n/ 2 * J 


1 


0 


$ n (x,z)e 1sz dz 


Similarly, 


(136) J "’ (x,s)5 if; j 


f0 *>,z)e 1sz dz 
38 


= cos 


Ai’xV 

-1 

l b J 

'JZTt{s~$ n ) 



and 


< i37 > V ; i' (x - s > ■ w f » <(X’^ isz * 

CTT 1 -oo . *\ 


- - r sin t- x < 


-i 


J27(s-B n ) , 


Fourier transforming (127) and (128) w.r.t. z gives 


2 

(138) ^ - Y 2 j3 n (x >s s) = 0 


^ 139 ^ (ax 2 - Y (K (X< ’ S) ° 


where y and y| are identical to (83) and (84). Also x> and x< 
have the same meaning as before. The general solutions to (138) 
and (139) are 


yX 

(1 40 ) $ n (x > ,s) 5 = $ n (x > ) = A(s)e" yX > + B(s)e > 

,(141) $ n (x <f s) = $ n ( x < ) = COsh y l x < + D ( s ) sinh Y 1 x < 


B(s) = 0 and D(s) = 0 from (133) and (130), respecti vely. Thus, 


(142) $ n (x >s s) = i+( X> ) * ;-( X> ) - A(s )e" yX 


(143) * n (x <# s) = £(*) + 3” (x < ) = C(s)cosh y,x . 


The steps leading to the final solution for $ n (x,z) are similar 
to those previously followed in the derivation of and s^^. 
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From (129) it follows that 


(144) ;+’( x =b ± )=0 


Using (132) together with (144) yields 


(145) -e^ y A(s)e~ Yb = Y-| C(s)sinh y-jb 


From (131) one obtains 


(146) ♦~(b ) = + $-( b ) 

" 72^(s - 6 n ) " ‘ 

(Eq. (136) has been used in obtaining (146)). One may re-write 
( 1 46) via (142) and (143) as 


(147) A(s)e' Yb -£(b.) = + c(s)cosh r,b - J*(b ) 

. " JRs - 6 n ) .. 1 " ‘ 


It is convenient to define 


(us) j + ( S ) s$*(b_> - ;*<b + ) 


Equation (147) becomes via (148) 


(149) A(s)e“ Yb - C(s)coshy,b = . j + ( s ) 

1 ./Ms - b„) 

Eliminating C(s) between (145) and (149) leads to 

(150) yA(s)e~ yb G~ ] (s) = — - j+( s ) 

J2 tT(s - B n ) 
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in which G(s) is identical to that in (94). One notes from (142) 
and (144) that A(s)e~Yb is given by 


(151) -yA(s)e" vb = ;;'(b + ) + $;'(b + ) = $-'(b t ) 


Utilizing (151) in (150) yields the Wiener-Hopf equation 

(152) G'’(s) ♦:'(b ) = J + (s) + 

" JM s - 6 n ) . . 


From the behavior of <j> n and 8<j> n /3x as |z| it is evident that 
(j>“'(b + ) is analytic for t < k? and J + (s) is analytic for t > - 
Min(k 2 ,k 4 ). Also G(s) is analytic for |tj<Min(k 2 ,k 4 ). Hence, one 
notes that Eq. (152) is analytic in the strip | T j < Min(k 2 ,k 4 ). 
Factorizing G(s) into G + (s)G (s) leads to 


(153) 


^'(b + )G: 1 (s) = J + (s)G + (s) + 


i (-1 ) n G + (s) 
- B p ) 


The L.H.S. of (153) is analytic in the L.H.P.; whereas, the R.H.S. 
of (153) is analytic in the U.H.P., except at s = $ n . For con- 
venience, (153) is re-written as 


1 K-l ) n G + (6_) + 

054) # n (b + )G_ (s) ±- 2 -= J + (s)G + (s) + 

Jm s - e n ) 

+ i(-l) n 


J2^(s - 3 n ) 


{G + (s)-G + (g n )> 


The L.H.S. of (154) is analytic in the L.H.P.; whereas the R.H.S. 
is analytic in the U.H.P. Thus, the two sides of (154) together 
define an entire function Q(s) which can be determined from the 
asymptotic behavior of (154) as |s| -»■ ». As before, G ± (s) <v 0(sl/ 2 ) 
as |s| -> ® in their respective regions of analyticity. The edge 
condition for this problem dictates that the electric field intensity 
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can become singular no faster than p' 1 ^ 2 in the neighborhood of 
the edge (at x = 0, z = 0) as p -*■ Q; since the transform of the 
electric field is proportional to the edge condition there- 
fore implies that (b + ) 0 (s _ V2) as | s | -► «, i n the strip. 

This is also deduced via techniques discussed in [9]. Thus, one 
concludes that Q(s) = 0 from the asymptotic conditions for large s. 5 
It follows that the L.H.S. and the R.H.S. of (154) are each equal 
to zero; hence. 


(155) 


5 = 


i(-i) n e + (e n )G_(s) 

^7(s - e n ) 


Equation (154) via (151) becomes 


(156) 


-YA(s)e" Yb = 


iH) n G + (e n )G_(s) 
- B n ) 



Also using (145) in (156) one gets 


(157) 



Y-j C(s)sinh, y-j b = 


i(-1) n G + (B n )G_(s) 
\/Ms - B n ) 


Incorporating the values of A(s) and C(s) from (156) and (157) into 
(142) and (143), respectively leads to 


(158) 

and 




-1 (-1 ) n G + ( 6 n )G_(S ) -y( X> -b) 

e 

y(s - B n ) 


(159) 3 n ( X< ,s) 


i(-D n e r G + (g n )Gjs) 

JTit{s - B n )Y 1 sinhy 1 b 


cosh Yt x 
1 < 



The formal solution to this problem is now complete since $ n (x>,z) 
and <f> n (x<,z) are respectively obtained from (158) and (159) by the 
inverse transformation. Thus, 
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(160) * n (x >s z) 


“ -i(-D n G + (e n )G_(s) -y(x > -b) ^ 

.00 - e n ) 6 6 ds 


and 


(161) * n ( x <’ z)= 57 


f i<-l) n E r G + ( 6 n )G_(s) 
. (s- 6 n ) Y,s1nh Y-jb 


—1 s z 

cosh Y-j x < e ds 


As before, the integrals in (160) and (161) may be evaluated via the 
Cauchy Residue theorem after closing the contours defined on the real 
axis by semi-circular contours at infinity. The semi-circular 
contours at infinity give vanishing contributions to the integral. 

The pertinent contours in the U.H.P. and the L.H.P. are indicated in 
Fig. 10. The fields for z > 0 are obtained by closing the contour 
of integration in the L.H.P. Whereas the fields for z < 0 are 
obtained by closing the contour of integration in the U.H.P. 


The field radiated into region ® is obtained from (160) 
via an asymptotic approximation to the integral. This radiation 
field arises from the branch cut contribution to (160). Introducing 
the polar transformations given by (108a, 108b) into (160), and 
deforming the transformed' contour of integration in the a plane 
to the steepest descent path provides one with the necessary form 
of the integral for a saddle point evaluation. The saddle point 
result which is valid for large kp yields an expression for the 
radiated field denoted by as 


(162) 

where 

(163) 


> rad (p,*) * F®4) ■ 

>/p” 


n (-1 ) n+1 G. (e n )G (-k sin *)e 

F®U) = + " ’ 


kb cos <f> + j 


(k sin ip + 3 n ) 


; n= 0 ,l , 2 ,- 


(Figure 8 a illustrates the coordinates p and <|> for the observation 
point. ) 

The TM 0 surface wave which is launched into regions (B) and® 
(for z < 0) due to the TM on waveguide mode incident at ® (from 
© ) may be obtained from the residue at the pole s = 8 of G_(s) 
in (160) for region ® (for z < 0 ), or from the residue at the 
pole s =3 of G_(s) in (161) for region ®. The residue cor- 
responding to the pole at s = 3 is obtained by closing the 
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contour of integration in the U.H.P. The result, for the (TNIq surface 
wave launched into region ® (for z < 0) is denoted by: where 


(164) = T^ w e “ a ^ x " b ^ e“ lSz ; x >.b,' z < 0 


in which T^ w is the surface wave launching coefficient given by 


(165) 


ySW_ 

n 


oG + (6)(6-6 n ) 




; n = 0, 1, 2, 


antT^Q was defined earlier in (113). One notes that G_(s) in (160) 
has been replaced by G(s)/G+(s) to arrive at (165). 

Finally, the fields reflected back into the waveguide region © 
are given by the sum of the residues corresponding to the poles at 
at s =-3 m of (y| sinh rib)" 1 in the integrand of (161 ); The 
residues at these poles are obtained by closing the contour of> >v i 
integration in the L.H.P. Before evaluating these residues, it 
is convenient to re-write (161) as 


-i(-l) n e r Mfg r G (s)cos 
( 166 ) $n(x< ,z)=— 


-°° J k d~ s2 sin J 


Jk 2 .-s Z x e" 1sz 
a — ds 


k^-s 2 b <S-y ; 


d 2 2 / ?2 -l 

kj-s sin>/kH-s b) cor- 
responding to fields in the region z >. 0 occur at s =-e m . Thus, 
the field scattered back into © which is denoted by is 
gi ven as . ' 


(167) = 


(-nVA) j o 


1 B m Z 

G (-3 )cos ^ x e 
- m b < 


d 

3? 




SI 


n^k^-s 2 b) 


£ 


(Qn+eJ 


m 


One may therefore express in (167) as 


<168) = jo R ™ cos ^ x < e1V 
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where denotes the reflection coefficient at © for an mth 
waveguide mode reflected back into region © ;due to an nth v. ; 
waveguide mode incident at © . R is thus given by 


(169) 


mn 


H)""" Er G + (6 n )G + ( 6m ) 


a„(B +6„)6„b 

mm n m 


where 


(170) 


A m = 


2, for m = 0 

1 , for m t 0 


(in going from (167) to (169), one again uses the relationship 
G_(-s) = G+(s) which is indicated in Appendix II). 

The generalized scattering matrix s -| 2 describes the fields 
produced in regions © and (jy by a TM on waveguide mode which is 
incident at © from the waveguide region © (i.e., from the right), 
The fields produced by the incident TM on mode consist of the 
radiation field in regions ® and and a TM 0 surface wave 
field which is launched into regions (§) and © (for z < 0). If 
one is interested in calculating the radiation field produced in 
region ® , then one defines s-j 2 to be Fjl(<j>) of (163); whereas, if- 
one is interested in calculating the TMq surface wave field which 
is launched into regions © (for z < 0) and (6), then s, 9 is 
defined to be T*w 0 f (165). Thus, u 


(171) 


’12 


= < 


T^ w , for calculating the surface wave field 

which_is launched into regions @ (for z<0) 
•and © - 

B 

F n ( 4> ) » for calculating the fields radiated into 
region ®. 


Clearly, s^ may be expressed as a 1 x M matrix where M -+• <*>. Thus, 


/ 1 75 \ e _ r T SW jSW*i*SW 

(172) s 12 = [T q ^ T 2 


r SW 


•] 1 xM ; for regions © (when z<0) 


and 


or 

(173) 


*12 


= F?U) f|(«) F®(t) ... fJJ(*) ...] )xH : 

for region ® 
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The generalized scattering matrix S 22 represents the modal 
coefficients of the fields reflected back into © due to a TM on 
waveguide mode incident at © from ©. The field reflected back 
into the waveguide is given by (168). Hence, it is clear the S22 
is an MxM, matrix (M -* °°) whose elements are R of (169). Thus, 


(174) s 22 = CR mn ] MxM 


^0 R 01 R 02 ** j 

R 10 R 11 R 12 “* 

• • • 




MxM 


Determination of the surface wave 
diffraction coefficient by the GSMT 


Referring to Fig. 9a, it is seen that the incident TM 0 surface \ 
wave field is scattered by junction © to produce a diffracted field 
in region ® ; one denotes its magnetic field by uS*, where uV is 
given via (122) as 


(175) 


ik 


s n u 


■JV 


p -AwS e 


ike 


vr 


; in region ® 


u 1 in (175) denotes the magnetic field of the incident surface wave 
evaluated at CD. The part of the incident surface wave field which 
is scattered by © into region © becomes incident at © , from 
which it is reflected; this reflected field is incident back at © 
where it undergoes further scattering into regions ®, © and ©, 
and so on, thereby giving rise to multiple interactions between 
junctions © and (f). The fields resulting from these multiple 
interactions may be expressed, in a Convergent Neumann series as 
done by Pace and Mittra[7]; however, an alternate procedure based 
on a self consistent method leads directly to the same result. 

The self consistent method will be used in this analysis; the 
use of this method is described with the aid of Fig. 11, Let u-]? 
represent the value of the magnetic field incident at © from (2) 
after taking all the multiple interactions into consideration. 
Similarly, let U 2 ] represent the value of the magnetic field 'I 
intensity incident at © from © after taking all the multiple 
interactions, into account. The total diffracted field in region © 
denoted by u^ is thus a superposition of the field u?., and u“, the 
field diffracted into © due to u -| 2 incident on ©. Thus, 


(176) u d = u d + u d 
where u d may be expressed as 
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rl R kp 

(177) u° = s 12 u 12 — = [F“]u 12 — ; in region ®. 

•I p J p 

R B 

s 12 = [ F n ] 1,1 t his case and F n is defined in (163). The expressions 
for u-|2 a nd U 2 i for 6 -► 0 (see Fig. 11) are given in terms of s 21 
and s 22 as 


(1/d) u 2 i s 2 i u + ^22^12 

(179) u 12 =r U 21 


where f is the generalized matrix for junction ©; it corresponds to 
the reflection coefficient at Q) and it is an identity matrix of 
infinite order I as mentioned earlier. Eliminating ui2 between (178) 
and (179) yields 


(I 


^22 ) u 21 ^°" 1 ^ 


21 


Hence it follows that 


(180) 


2.1 


= (! - s 22 ) 


-1. 


5 21 


From (179) and (180) u^ 2 becomes 
(181) u 12 = (I - s 22 )" 1 s 21 u 1 - 


Incorporating (181) in (177) yields 


(182) ujj - s 12 (I - s 22 ) 1 s 21 u 1 

Finally, combining (175) and (182) according to (176) yields 
the expression for the diffracted field u d as 


(183) 


U = [s-j-j + S 15 (I - S99) s ?1 ju 


ikp 


’ 12 ' 


22 ' 


21 
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REGION (A) 


U'=TM 0 SURFACE 
WAVE 

FIELD AT © 


m 


JUNCTION 


(b 


— s — 

REGION © 


U 2 , U 2 | W 

l 

U I2 ^ U 2I ^ ' ■ [ ‘ 

JUNCTION (2) 


Fig. 11— Multiple interactions between junctions 0 and Q . ) 
for the calculation of the diffracted field. 


Eq. (183) may be written as 

a ,• Jkp 

(184) u d = D 2 u 1 p 

•ip 

where D 2 is the diffraction coefficient for the canonical problem of 
the recessed, truncated dielectric cover which is flush mounted in 
a perfectly-conducting surface (see Fig. 2). Thus, the formal 
expression for D 2 in terms of the scattering matrices of junctions 
0 and 0 is obtained via (183) and (184). 


(185) D 2 - S-ji + s 1 2 ( i ” Sgg) Sgi 

A B 

Since sn = F (<*>), and s 12 = C ) ], 085) becomes 

(186) D 2 = F A (<j>) + [f|J(,|>)](I - s 22 r 1 s 21 

where s 22 and s 2 -j are given by (174) and (123), respectively. 
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Determination of the surface wave 
reflection coefficient by the GSMT 


The expression for the surface wave reflection coefficient for 

the canonical problem illustrated in Fig. 4 is derived in a manner 

analogous to that for the diffraction coefficient discussed previously. 

Figure 12 illustrates the self consistent procedure for determining 

the reflected TM 0 surface wave field for the canonical problem. 

The reflected TM 0 surface wave field is denoted by u r . The incident 

field at © is denoted by u 1 as before. Also, u-j 2 and 1*2 1 used 

here have the same meaning as in the previous discussion for the 

diffraction. coefficient. The field u r is a superposition of the 

field expressible in terms of the scattering matrix sn at 0 

denoted by u-T, and the field contribution due to multiple interactions 

denoted by u^. Thus 
m 

(187) u r = u\ + u r m 


where 

(188) U 1 f sw 


and 


(189) ^ = s 12 u 12 f sw 


The quantity f in (188) and (189) is defined by 


(190) 


-a(x-b) e -i j n region (§) (for z<0) 




cos V - 1 ez 

” 9 


cos k x b 


in region ®. 


The quantities a, B and k x in the expression for f sw are the same as 
in (62) and (63). Following the procedure for the diffracted field, 
one may write 
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REGION (A) ( FOR Z<0) 
OR REGION (B) 


REGION © 




'm 


U 2I U 2I' 

■ U 12 Tu 21- 


JUNCTION 


d> 


JUNCTION 


<t 


Fig. 12-- Multiple interaction between junctions ® and ® 
for calculating the reflected surface wave field. 


(191) 


J 12 


= (I 


’22 


) 1 s 21 u 1 


so that (189) becomes 

092) ujj = s 12 (I - s 22 ) 's^u’f^. 


and u r becomes 

(193) u = ( s 1 1 + s-] 2 ( I - s 22 ) s 2 i } u"f sw 


The above expression may be written as 
(194) u r = uV sw 

where R 2 is the TM 0 surface wave reflection coefficient for the 
canonical problem of the recessed, truncated dielectric slab which 
is flush mounted in a perfectly conducting surface (see Fig. 2). 
f *2 is given via (193) and (194) as 
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(195) Rj E R S " + [T= W ] (I - s 22 ) _1 s 2 , 

because s-j -| = R sw and s-j 2 = f° r the surface wave reflection 

coefficient calculation. R sw and T^ w are defined in (115) and 
(165), respectively. The quantities s ?? and s ?1 are given by (174) 
and (123), respectively. 



Fig. 13-- Geometry for the variational solution. 


Solution Based on the Variational Technit 


In this section, a variational solution to the canonical problem 
is briefly indicated. The geometry for the problem is indicated in 
Fig. 13 above; it is identical to Fig. 2, except for a shift in the 
origin of the x-z coordinate system. The region defined by (x > -b 
and z < 0) is denoted as region I; whereas, the region for (x >~D, 
z > 0) is denoted as region II. The excitation for this problem is 
again a TM 0 surface wave incident at (x = 0, z = 0) from region I 
(i.e., from the left). The incident and the total magnetic fields 
are denoted by and $. Here, 0 = in region I; whereas, $ = 
in region II. 

f 

-ox igz 

e ; for x > 0 and z < 0 

(196) *1 

cos k (x+b) 

— - — £-7- e e ; for -b < x < 0, and z < 0 

cos kb — — 

X 

\ 


where a, b and k x are identical to those in (62) and (63). One may 
express the fields and in terms of special Green's functions 
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Gj and GII for regions I and II, respectively, as.., 


••• ' . . yv- f “ 

(197) <fj(x,z) = ${(x,z) - G.(x,z;x' ,0) -jp 

Jq 


dx' + ($;{(x,z))*, 
• =n 1 


z'=0 


for region I, and 


f°° 3 r 

(198) *jj(x,z) = G n (x,z;x* ,0) —f 

o 


z'=0 


dx* , 


for region II.. The term ($}(x,z))* in (197) refers to the complex 
conjugate of $|(x,z) and is present as a result of the boundary 
condition given in (200) so that in effect there is a perfectly- 
conducting wall at z = 0, x > -b. The special Green's functions 
Gj and Gj i satisfy the following differential equations and boundary 
conditions. 


(199a) + + k 2 ) Gj(x,z;x‘ ,z' ) = -6(x-x' 

\ 9 X 9 Z J 


)6(z-z' ); 
for x > 0, and z < 0. 


(199b) 


{ 


9 2 9 2 2 

9/ 9? d 


lj Gj(x,z;x' ,z' ) = 


0 ; 


for -b < x < 0, and z < 0 


( 200 ) 


9G 


9Z 


= 0 


z=0 


( 201 ) 


3G 


9X 


= 0 


x=-b 


(202) 0,1 


(203) 


gGj 

"ax" 


= 1 9G I 

x=0+ e r 3x 


x=0- 
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(204) G, satisfies the radiation condition as | (x 2 +z^)^ 2 |-*>° 

* «j ^ 

for an e” U)L time dependence. 

The quantities, k, ; kj and e r appearing above are the same as those 
in (62) and (63). The 6 in (199a) symbolizes the Dirac delta 
function. 


(205) 


(206) 


( 


a 2 a 2 2 

^T +k 
ax az 


) ^1 j(x,z,x , 


V) = 


-6(x-x' )s(z-z' ) ; 

for x > 0, and z > 0. 


aG 

az 


II 


= o 


'z=0 


(207) : 


aG 

ax 


II 


'= 0 


r x=0 


(208) G tt satisfies the radiation condition as | (x 2 +z 2 )^ 2 |-*>° 

for an e" time dependence. 

The derivations for the special Green's functions Gj and Gj j 
will not be presented here for the sake of brevity; only the final 
forms of Gj and Gjj are given. 

(209) Gj = - 6 G x (x,x') G z (z,z')dx xo 

c x 

where 


(210) G x (x,x') 


and 

-i Ik z cos /a, z 

p 'I z < V z > 

(211) G ( z , z * ) = = , for z < 0 

-U x z 
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Here, x< refers to values of x for which x < x 1 , and the quantities 
x >» z > and z< are defined in the same manner. 


( 212 ) 


x xo + x z - k 


and 

(213) 

with 

(214) 


- r £ rfxo cos fx1 b + i J x x1 sin J X xl b 

£ rj*7o cos fTl b * *A7l s1n J^l b 


x xl + x z = k d 


The contour C x in (209) encloses the singularities of G x (x,x' ) in 

5 

i u | x-x ' | + iu(x+x') i'Jk 2 -u^z > 


the complex A plane. Also, 

AU 


(215) 


’II " 2wi 


du 




cos 


J^ z< 


The Green's functions Gj and Gn given above have been constructed 
via the techniques indicated in [12,13]. The integral in (209) 
may be re-written as a superposition of a pole contribution 
and a branch cut contribution associated with the singularities of 
G x ; the pole contribution may be expressed in the form R* (o](x,z))* 
in which the complex constant R' is unknown at present; whereas, the 
branch cut contribution has the form (after employing a change of 
variables u 2 = a X q) 


(216) 


2*i 


dx' 


iu|x-x' | + ^ e iu(x+x‘ ) 


J o 



> - U2 >Jz’=0 


9Z • 


z'=0 


54 



Hence, one may rewrite (197) and (198) as 


(217) 


(x,z ) = (x ,z ) + (1-R‘ )($](x,z))* - 



du 


■*0 


dx'S 


f \ 


r f 

^ e 1 u i x-x ' 1 + r e iu(x+x' ) 

* • 

, 'N 
(e 

, 7k 2 -u 2 



? 2 9<j> t 

• cos yj k -u z > jp- 


J z'=0 


and 

(218) 


»ii (x > z) 


,oo fOO 


du 


dx' i 


jo 


; iu|x-x* | + e iu(x+x' ) 

l ” 


WP3 


cos 


’ < 


az’ 


J z'=0 


The fields and $n must be continuous across the aperture defined 
by z = 0 and x > 0. Also 9$t/9z' Li n and 9 $ TT / 9 z' |. n must be 
continuous. Le^ 2-0 u 2 " u 


(219) 


34>j 

9 Z 1 


9$ 


II 


z'=0 


9 Z 


z'=0 


K(x') 


Equating $j(x,0) to $n(x,0) and utilizing (219) one obtains the 
following integral equation 


( 220 ) 


(1 + R 3 )$j(x,0) = 


* 

r(x,x' )K(x' )dx' 

0 


where 
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( 221 ) 

and 

( 222 ) 



- R 1 


r(x,x' ) 


1 

2 Tri 


; e i|x-x'[ t ( 1 4 )e 1 u(x*x l ) 

J7T7 


du 


The quantities R 3 and K(x') in (220) are unknowns. It is clear 
from (217), that R 3 as defined by (221) represents the surface 
wave reflection coefficient which is of interest. Multiplying 
both sides of ( 220 ) by K(x) and integrating on x gives 


(223) 


(1 + r 3 ) 


* • 

K(x)$!(x, 0 )dx = 
JO 1 


CO #00 

K(x)r(x,x' )K(x' )dx' dx 
0J 0 


One may represent K(x) by a sum of discrete surface wave modes 
(arising from the discrete spectrum of the integral operator in 
( 220 )) and an appropriate integral (pertaining to the continuous 
spectrum contribution); hence, 

( 224 ) K(x) = |jr [♦ 1 I (x,2 , )+R3*{*(x,z , )]i I ,. , 0 +» c (x,6) 

= 1 s( 1 -R 3 )tJ(x, 0 ) + * c (x, 0 ) . 

The present analysis is restricted so that obly the dominant surface 
wave mode exists. The term $c( x >°) denotes the contribution to ' 
K(x) arising from the continuous spectrum. The surface wave modes 
are orthogonal to the continuous spectrum; thus (224) when multiplied 
by $|(x, 0 ) and integrated with respect to x gives 

00 

(225a) (1-R 3 )=ff| K(x)$j(x, 0 )dx A 


in which 


(225b) 


00 


|V T (x, 0 )] 2 dx = 
J 0 1 


(e“ aX ) Z dx = 
j 0 


1 

27 


has been used. From (223) and (225a) one obtains the following 
stationary form for (l+R ^)/ (1 -R3) 
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f CO #co 


(226a) 



K(x)r(x,x' )K(x' )dx' dx 


i e ' 0*1 

0 

2 “T 

K{x) $!(x,0)dx 

Jo 1 J 

2 


To solve for R3 in (226a), it is customary to assume that K(x) can 
be approximated by the incident and 1 reflected dominant surface wave 
modes, i.e., 

K(x) £ ig(l-R 3 )$ 1 I (x,0) . 

When the above approximation is chosen as the trial field, (226a) 
becomes 


(226b) 


1 + R, r m 


0 


fCO 

dx j dx' (x,0)r(x,x' )<J>j(x' ,0) 


The integral on the R.H.S. of (226b) can be simplified by performing 
the integrations on x and x'. 


(227) 


r oo r co 

dx 
0 


J 


dx' e“ aX r(x,x' )e~ aX ' =^P.V. 


1 + R 


T 


L a fc + U* (a - i U ) 2 J 


du 


y . j 


The symbol P.V. in front of the integral on the R.H.S. of (227) 
indicates the principal value of the integral which is evaluated 
numerically- to calculate R3. The evaluation of the diffracted field 
is discussed next. The branch cut contribution to 



gives rise to the diffracted field in region I denoted by $j. Thus 


(228) 
*1 = 


1 1 

2tt1 


du 


J-c 


dx' 


+ ft e 1 * u ( x+x ' ) V -t J k2-u 2 z 

1 y~j ‘ 


K(x') 
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The branch cut contribution to the field corresponds to the 
diffracted field in region II; the diffracted field in region II is" 
denoted by <^j. Thus, S. 


(229) 

,d _ T 
$ II ' 2m* 


du 


0 dx ' ( ! 


/ iulx-x' | iu(x+x' )\ . I. 2 2, 

• I e 1 ■ e l__j \ e Wk -u z K(x , , 




2 2 
k - u^ 


Employing K(x‘ ) £ i3(l-R3)${(x',0) in (228) and (229), and inte- 
grating on x' yields 


(230) 
*1 “ ■ 


1 

2r\ 


du 


e 1ux e^ 

r 2 2 

>jk. - u 


,2 2 , 

k -u z 


i 3 ( 1 - R ) 


'ja + i U + a - i U | J 


and 


(231) 



1 

2-rri 


• L2 2 


r > 


r°° ux ilk -u z 

Hu e e ^ 

i 3(1-R) - 

2a ! 

r 

au 

2 2 1 


Jk 2 - u 2 

- 

a + U 1 

. 


One may evaluate (230) and (231) asymptotically for large kp via the 
saddle point method. As before, one introduces the polar trans- 
formation 

(232) u = k sin E, , J k^ - u^ = k cos £ 
and 

(233) x = p cos <f>, z = p sin 4 > . - . 


The coordinates (p,4>) at the observation point are indicated in 
Fig. 13. 
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(234) 

id 


*1 = - 


1 

2ri 


i*k P sin(?-4>) 


8(1-R 3 ) 


f 


+ i 


+ R(k sin g) 

k sine a - ik sind 


de 


The contour in the complex e plane is deformed into the steepest 
descent path, and the integral is then evaluated via the saddle 
point method to obtain 


(235a) 


i en- | V A o 

J 2iirk 


1 


a + ik COS <f> 


+ 




R(k cos <t>) l 
a^nr~co?i’ f 



-ir/2<(f)<0 


Similarly, the saddle point result for $ 


d 

II 


is 


(235b) 



„ m - b.i) 

j2Uk 


:a 


?? 2 
a + k COS (j> 



0 < (p < tt/2 


The constant Aq is introduced so that $t( 3>=0) = <j>n( =0); there is 
no other analytical justification for its presence. The quantity 
R3 in (235a) and (235b) is evaluated from (226b). Clearly, the 
quantities which multiply the e^P/Tp" factor in (235a) and (235b) 
correspond to the surface wave diffraction coefficient for the 
regions -tt/2 < - 4> < 0, and 0 < <f> < w/2, respectively; this surface 
wave diffraction coefficient which is based on the variational 

solution is denoted by D 0 . 
r -j 


(236) 

°3 E i 


sj2Uk 

ib(i-r 3 ) 
\j2UF 


1 + R(k cos 4>) 

a + ik COS <)) a - ik COS <j> 


2 2 2 
a + k COS 


-tt/2 < < 0 


0 < <p < tt/2 
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CHAPTER IV 

RESULTS AND DISCUSSION 


The TM 0 surface wave reflection and diffraction coefficients 
associated with the termination of the dielectric cover in Fig. 2 ■*'[ 
have been obtained in this report; these coefficients allow one to 
employ the GTD method to calculate the fields radiated by a con- 
figuration such as that shown in Fig. 1. The GTD method is chosen ' 1 
because r 

a) it is simple to use, and yields accurate results 
in high-frequency problems; 

b) it provides some physical insight into the radiation 
and scattering mechanisms involved; and 

c) it can be used to treat problems for which exact 

analytical solutions are unavailable. 

Furthermore, the GTD method has been employed successfully in the 
past to calculate the radiation from a dielectric covered slot in 
a finite ground plane in which the dielectric completely covers 
the ground plane [14], thereby providing an added confidence in the 
method. 

In the GTD method, the total field at P is obtained via a 
superposition of the geometrical optics field along the ray s 0 , 
and the diffracted field at P which arrives via the ray paths si 
and S 2 , as illustrated in Fig. 1. As in Chapters II and III the 
y-directed magnetic field is denoted by u. Let the geometrical 
optics field be denoted by u g o.» and the diffracted field from 
the terminations A and B by ug and ug, respectively. u g>0 cor- 
responds to the direct radiation from the slot to P given 6y 


(237) 




The slot also excites surface waves which diffract from the termi- 
nations A and B. The surface wave excited by the slot is denoted 
by u sw , where 


(238) 


% = K o e 


±1 3Z 



for z 


> 

< 


0 
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in which a and g are the surface wave attenuation and propagation 
constants discussed previously in Chapters II and III. The 
quantities C 0 and Ko are deduced from the solution for the fields 
of a narrow slot radiating through a dielectric covered plane of 
infinite extent, as shown in Fig. 14. The quantities C 0 and Kq 
will not be derived here since they may be obtained via techniques 
similar to those outlined in [8] for an electric line source ex- 
citation of a grounded dielectric slab. Only the final results 
for C and K are presented below. 


(239) 


C. = 


i we M„ cos <p e 
oo 


■'i 


and 


J 2 ’ k {£ -sin% sin kje^sin^b + i e r cos $ cos k 


e r -sin*<|,b] 


(240) 


“ ME o k x a M o 

agO+e^bJsin k x b + k x g(ab+e r )cos k x b 


2 2 2 

The quantities k x and b correspond to (e r k - g ), and the dielectric 
slab thickness (in meters), respectively. The quantity Mg in (239) 
and (240) refers to the strength of the equivalent magnetic line 
source associated with the electric field in the aperture of the narrow 
axial slot. The total field radiated at P by the configuration in Fig. 
1 is obtained by superposing the geometrical optics contribution in 
(237) and the diffracted fields u^ and ug, respectively, where 

K e m ° iks 1 

(24U) ufo). ° i2 

1-R e 0 JS 1 


and 

(241b) 


’Wq 

U 8< P > * 

1-R e c 



for the case when the narrow axial slot is equidistant from the 
terminations at A and B (the distance being £ 0 ); the corresponding 
expression for the case when the slot is not equidistant from A 
and B is not presented here, but it can be derived in a straight- 
forward manner. In Eqs. (241a) and (241b) the contributions from 
all surface waves multiply reflected between terminations A and B 
have been included. The quantities R and D in (241a, b) refer to the 

TM„ surface wave reflection and diffraction coefficients associated 
o 
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with the terminations at A and B, respectively. Thus, the total 
field, u(P) at P is 

(242) u(P) = Vo> ♦ u d ♦ u B d 

where Ug. 0> , and ,Ug are given in (237) and (241a, b).. Although 
examples of this type are two dimensional , the patterns obtained are 
identical with those of a related three-dimensional problem, in a 
plane of symmetry perpendicular to the axis of the slot. 



Fig. 14 ---Narrow axial slot in an infinite ground plane 
covered with a dielectric slab of infinite 
extent in the z direction. 


As indicated in the introduction, the quantities J R and D in 
(241a, b) are derived by three different methods, namely: 

1 . M' ' • 

(1) An exact solution to the canonical problem of Fig. 2 
based on the GSMT; the details of this solution are given 
in Chapter III. The expressions for R and D obtained 

by this method are denoted by R 2 and D 2 , respectively. 

R 2 is given in (195) and D 2 is given in (186); 

(2) An approximate solution to the canonical problem of 
Fig. 2 obtained by the variational method; the details 
of this solution are given in Chapter III. The ex- 
pressions for R and D obtained by this method are 
denoted by R 3 and D 3 , respectively. R 3 is given in 
(226b) and is given in (236); 
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(3) an exact solution of the problem pertaining to the 
approximate physical configuration of Fig. 3 in which 
a reactive surface approximates the dielectric cover 
of the exact configuration of Fig. 2; the details of - 
this solution are given in Chapter II, and the ex- 
- pressions for R and D obtained by this method are 
denoted by Ri and Di , respectively. R-| is given 
in (53) and D-j is given in (60). 

The numerical results for the TM 0 surface wave reflection 
and diffraction coefficients R and D, respectively, which are 
calculated by the three different methods are presented in this 
chapter. A comparison of the numerical results based on the 
different methods is made to estimate the accuracy of the 
approximate solutions. The numerical results for the TMq 
surface wave reflection' coefficient are given first; these are 
followed by calculated curves for the TM surface wave diffraction 
coefficient. 

We wish to compare the reflection and diffraction coefficients of 
the terminated surface reactance model with those of terminated dielec- 
tric slab structure of Fig. 2. This is made possible by the relation- 
ship between X s and b' for a given e*., see Eqs. (2) and (5). Thus, 
in the following R-j and D-j are given both as functions of X s and 
b (for a fixed e r )- Calculated values of the surface wave reflection 
coefficient R] are listed in Table I as a function of the normal 
surface reactance. The thickness of the grounded dielectric cover 
with the same normalized surface reactance is listed in column two. 

The values of |R-|| given in Table I are in agreement with those ob- 
tained from Kay's solution [3] to the problem in Fig. 3; however, 
the phase of the reflection coefficient in Kay's solution is not in 
a form tractable .for numerical computations, and therefore it cannot 
be compared with the phase of R-j presented in Table I. As mentioned 
previously it appears that the only reference which contains a solution 
from which phase.calculations are tractable is that of Weinstein [5]. 
However, the results in [5] have an error in the final form; the error 
is corrected fairly easily and the numerical values for the reflection 
coefficient based on the corrected results are included for com- 
parison. Both the magnitude and the phase of the reflection coefficient 
calculated from [5] are also given in Table I; it is seen from there 
that the values calculated from the corrected results in [5] agree 
well with the values of Ri. The magnitudes are in complete agree- 
ment; only the phase of the reflection coefficients calculated by these 
two methods differ slightly. Thus, we have a valuable mutual check 
between the solution given in Chapter II, and the one in [5] obtained 
by a different factorization method. 
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The value of the reflection coefficient R 2 obtained by the 
formally exact GSMT solution for the configuration of Fig. 2 are 
given in Table II. These values are also presented graphically 
in Fig. 15. The GSMT result for R 2 involves infinite order 
matrices which are truncated to finite orders for numerical cal- 
culations; it is found that accurate computations are possible 
with only moderate size matrices because the results converge 
rapidly with an increase in matrix size. In these numerical 
calculations, the largest matrix inverted is 5x5; in fact, for 
all practical purposes, most of the results obtained here are 
found to converge with only a 3 x 3 matrix. In some cases, an 
even smaller matrix is adequate. The convergence of R 2 as a 
function of the matrix size is shown in Table III for several 
cases. 


TABLE II 

Reflection Coefficient R 2 (for the problem in Fig. 2) 
Based on the GSMT Solution 
e = 2.49, f = 8.9 GHz, 
free space wavelength x Q = 0.0337 meters 


Thickness of the 
dielectric cover, 
b (meters) 

Magnitude 
of R 2 

Phase of R£ 
(degrees) c 

.001 

.01545 

-8.8990 

.002 

.06955 

-15.878 

.003 

.18441 

-20.758 

.004 

. 38495 

-20.978 

.005 

.62913 

-15.309 

.006 

.80077 

- 8.643 

.007 

. 88505 

- 4.674 

.008 

..9260 

- 2.734 

.009 

.9483 

- 1.747 


( Note : A plot of the values in Table II is 

shown in Fig. 15). 


Values of the surface wave reflection coefficient R 3 calculated 
from an approximate solution by the variational method for the exact 
canonical configuration of Fig. 2 are listed in Table IV. It is 
seen from Tables II and IV that the values of | R 3 1 calculated by the 
variational method are consistently lower than tne corresponding 
values of ] R 2 | which are based on the exact solution. Similarly, 
Hwang [17] has found that a variational solution for the problem 
illustrated in Fig. 3 also gives numerical values for the magnitude 
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(b) 

Fig. 15— Magnitude and phase of based on GSMT solution. 
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TABLE III 


Convergence of R2 as a Function of the Matrix Size 
e r = 2.49, f = 8.9 GHz, free space wavelength a q = 0.0337 meters 


Thickness of 
the dielectric 
cover b 
(meters) 

Order of the 
largest matrix 
inverted 
N 

Magnitude of R£ 

Phase of R? 
(degrees) 


1 

.015362 

-8.5226 


2 

.015430 

-8.8028 

.001 

3 

.015444 

-8.8678 


4 

.015449 

-8.8946 


5 

.015450 

-8.8990 


1 

.346 8 

-23.926 


2 

.3777 

-21.478 

.004 

3 

.3824 

-21.138 


4 

.3841 

-21.028 


5 

.38495 

-20.978 


1 

.83465 

-15.899 


2 

.87939 

- 6.188 

.007 

3 

.88319 

- 5.148 


4 

.88445 

- 4.821 


5 

.88505 

- 4.674 


TABLE IV 


Reflection Coefficient R3 (for the problem in Fig. 2) 
based on the variational solution 
e r = 2.49, f = 8.9 GHz, 
free space wavelength \ Q = 0.0337 meters 


Thickness of the 
dielectric cover, 
b (meters) 

Magnitude 
of R 3 


.001 

.0072 

-19.02 

.002 

.0349 

-30.80 

.004 

.2061 

-37.37 

.006 

.4320 

-23.93 

.007 

.4812 

-18.76 

.009 

r- 

.5105 

■ \ 

-13.76 
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of the reflection coefficient which are consistently lower than the 
corresponding values of |R-]|; a graphical comparison of the values 
of 1 R-j | with the variational calculations of Hwang is made in 
Fig. 16. 



X s =O.1 128 X s = 0.3581 X s = 0.6245 X s =0.9047 X s = 1.035 


Fig. 16— A comparison of the variational and exact calculations 
for the reflection coefficient pertaining to Fig. 3. 
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The numerical results for the TM 0 surface wave diffraction 
coefficients are given next. The magnitude of Di for the con- 
figuration in Fig. 3 is plotted in Fig. 17a as a function of the 
aspect angle <f>, and for different thicknesses of the dielectric 
cover (or equivalent surface reactance). Corresponding numerical 
results for the phase of Dl are shown in Fig. 17b. The values of 
| Di ] agree with the corresponding values which are calculated from 
Kay's result [3] for the surface wave diffraction coefficient; how- 
ever, the expression for the phase of the diffraction coefficient 
in [3] is in a form intractable for numerical calculation and hence 
cannot be compared with the values of the phase of D] (shown in 
Fig. 17b). 

The numerical values of the surface wave diffraction coefficient, 
p2 calculated by the GSMT for the configuration of Fig. 2 are shown 
in Figs. 18 through 22. These values are normalized such that 0 dB 
corresponds to the maximum value of the quantity 20 logio ^id<|D2| 
for the thinnest dielectric cover considered; namely for the 
dielectric cover whose thickness, b = 0.001 m. Also, the 0 dB 
level in Fig. 17a is identical to that in Figs. 18-22. Curves for 
the phase of D2 as a function of aspect are plotted in Figs. 23 through 
26. The formal result for D2 which is based on the GSMT solution 
also involves infinite order matrices which must be truncated to 
finite orders for numerical calculations. Figures 20, 22 and 25 
show the convergence of the calculated values of D2 as a function of 
the size of the matrix being inverted; indeed, it is seen from these 
results that the convergence is very rapid. The R2 and D2 calculations 
require an increasing matrix size whenever the thickness and/or the 
permittivity of the dielectric cover increases. 

It is noted that the calculation of the field in the neighborhood 
of the edge (x = 0, z = 0 in Fig. 2) requires an increasing number of 
modes, i.e., an. increasing matrix size; in fact the proper edge be- 
havior can result only if one includes an infinite number of modes 
(i.e., if one uses infinite order matrices) as indicated in [10]. 
Fortunately, in this problem, one is interested only in the far-zone 
radiation field, and the reflected surface wave; both of which are 
evaluated sufficiently far from the edge. 

Figures 27, 28 and 29 give a comparison of the numerical values 
of I D] | , ] D2 I and | D3 J for a fixed thickness of the dielectric 
cover. It is seen from Figs. 27, 28, and 29 that the values of 
| Di | and |Dj| agree well with the values of | Do | i n the forward 
aspects defined by 0 < <j> < tt/2, which is also the region of higher 
field levels; the differences between | D2 | and j Dq [ are approxi- 
mately 6 to 10 dB in the backward region of space defined by 
-tt/2 < 4 > < 0. It is seen that the agreement between | D3 1 and | D2 1 
improves as the thickness of the dielectric cover increases; this 
occurs because the diffracted field along the aperture (z = 0, 
x > 0 in Fig. 13) becomes weaker with an increase in the thickness 
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of the dielectric cover, thereby making the trial field in the 
variational solution correspondingly more accurate for pattern 
calculations. The strength of the diffracted field becomes 
weaker with an increase in the thickness of the dielectric cover 
because the surface wave is now more tightly bound to the 
dielectric-air interface, thereby increasing the magnitude of 
the reflection coefficient; this in turn results in a corresponding 
decrease in the energy diffracted by the termination. 

It is seen that the diffraction coefficient D 3 is more 
accurate than the diffraction coefficient Di when the latter is used 
to approximate D 2 . The surface wave diffraction coefficient is 
proportional to the radiation pattern function, and hence Figs. 18- 
22 describe the radiation patterns associated with the diffraction 
of a surface wave by the termination of the dielectric cover. The 
nature of the radiation pattern is intimately connected with the 
value of the reflection coefficient of the surface wave reflected 
by the termination; this may be explained as follows. The radiation 
pattern may be considered as being produced by an equivalent magnetic 
current distribution located at x = 0 and for z <. 0. This current 
distribution is well -approximated by the incident and reflected 
surface wave fields at x = 0 and z <_ 0. One may therefore consider 
the equivalent current to be proportional to eisz + r e~i &z for 
z <. 0 where R is the surface wave reflection coefficient. The 
radiation pattern in the forward region is primarily produced by 
the termination of the incident surface wave at z = 0 ; whereas, 
the radiation pattern in the back region is primarily produced by 
the excitation of the reflected surface wave at z = 0. When R«l, 
the radiation pattern is expected to be much higher in the forward 
region than in the back region; whereas, when R-*l, the pattern 
levels are expected to be nearly equal in the two regions. This 
viewpoint is supported by the curves in Figs. 15 and 18-22. From 
the values in Tables I, II and IV, it is apparent that the surface 
wave reflection coefficients R 3 and R-j are not in close agreement 
with the values of R 2 . The variational solution is expected to be 
somewhat inaccurate because the present trial field cannot properly 
account for the significant diffraction effects, especially in the 
vicinity of the truncation of the dielectric slab; one notes that the 
present choice of the trial field in the aperture is equivalent to 
the commonly used Kirchhoff approximation (this approximation is 
also used by Angulo [ 6 ] in his variational solution for the problem 
in Fig. 4). One would expect the variational solution to be 
improved by incorporating a correction term to the surface wave 
fields in the aperture so as to provide some information on the 
fields scattered by the edge (x=0, z=0 in Fig. 2) along the aperture. 
Although the nature of this correction term is not known apriori, it 
could be based on the high frequency diffracted field in the aperture 
together with a modification which would provide the proper behavior 
in the vicinity of the edge. 
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Fig. 19 — Magnitude of the diffraction coefficient D. 
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Fig. 20-- Magnitude of the diffraction coefficient D, 



Fig. 21 — Magnitude of the diffraction coefficient D. 
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Fig. 22 — Magnitude of the diffraction coefficient D, 







Fig. 23— Phase of the diffraction coefficient D. 
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Fig. 24-- Phase of the diffraction coefficient D, 



Fig„ 25 — Phase of the diffraction coefficient D. 



Fig. 26 --Phase of the diffraction coefficient D, 
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Fig. 27 — A comparison of the diffraction coefficients 
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Fig. 28 — A comparison of the diffraction coefficients 



Fig. 29— A comparison of the diffraction coefficients 



The inaccuracy in the values of D-j and R] for the surface 
reactance approximation to the grounded dielectric cover of 
Fig. 2, appears to.be due largely to the fact that it does not 
properly account for the effect of the conducting step at z=0 
and 0 < x < b in Fig. 2. As mentioned in the introduction, the 
surface reactance approximation indicated in Fig. 3 also serves 
to approximately describe the diffraction of a surface wave by 
the complementary structure of Fig. 4. It is conjectured that 
the terminated surface reactance model, of Fig. 3 approximates 
the configuration of Fig. 4 better than it approximates the 
configuration of Fig. 2, because the latter structure has an 
additional discontinuity in the form of the step in the ground 
plane; as mentioned previously, this ; is evident by noting that 
the former structure reduces to the smooth ground plane as 
whereas the latter structure has a residual step present. The 
importance of this residual step on the radiation properties of 
the junction may be demonstrated as follows. If one maintains 
the effective surface reactance of the grounded dielectric 
cover constant while* decreasing the thickness of the slab for 
the configuration of Fig. 2, one tends to reduce the effect of 
the step without changing the discontinuity in surface reactance. 
The time average power ratio of the fields outside the dielectric 
to the fields inside the dielectric increases rapidly as the step 
size is decreased, while maintaining the surface reactance fixed 
thereby substantiating the argument that the effect of the step 
is reduced as the step size is decreased. Under these conditions 
one would expect the results for the configuration in Fig. 2 to 
approach those corresponding to the surface reactance configuration 
of Fig. 3. Indeed, this is found to be the case as shown in 
Figs. 30(a) and 31(a), where it is seen that | R 2 1 approaches |R-|| 
for a given surface reactance, X s , as the step size decreases. 

One notes that the relative permittivity, z? of the dielectric 
cover increases with diminishing step size in order to keep the 
equivalent surface reactance a constant. The phase of R 2 also 
approaches the phase of R] until a certain minimum step size b 0 
is reached, see Figs. 30 and 31; however, for any further decrease 
in the step size, the phase of R 2 tends to deviate from the phase 
of R-|. This anomalous behavior for the phase of R 2 is not in- 
dicated in Figs. 30(a) and 31(a); it may result from several 
factors which are considered next. The phase anomaly occurs 
whenever the permittivity (e r ) of the dielectric slab increases 
most rapidly with decrease in the thickness, b of the dielectric 
slab. It is seen that z r increases extremely rapidly with only 
a slight change in b for 0 < b < b 0 . One reason for the phase 
anomaly may be the possible existence of undetected numerical 
inaccuracies in the computer programs for evaluating the Wiener- 
Hopf factors of G(s) of (94) whenever e r becomes very large (the 
factors of G(s) and their numerical evaluation is discussed in-. 
Appendix II). Another factor contributing to the phase anomaly 
may arise from the following possibility. It is found that as 
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Fig. 30 — Behavior of the reflection coefficient FL 
function of b, for X„ = 0.1125. L 
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the height b of the step is decreased, the electrical thickness of 
the dielectric slab decreases only slightly and then actually 
increases; this increase in the electrical thickness becomes most 
rapid where e r begins to increase most rapidly. It is possible 
that when the electrical thickness of the step increases rapidly, 
one may not be able to isolate the junction effects into those 
which are assumed to arise independently from the discontinuity 
in surface reactance and from the presence of the step. 

In summary, it is seen that approximations based on the 
surface reactance model and the variational solution, where the 
trial function consists of only the surface wave fields, do not 
yield accurate values of the reflection coefficient for the 
canonical problem shown in Fig. 2. On the other hand, these 
approximations do predict the radiation pattern in the forward 
region with reasonable accuracy, but they fail to yield accurate 
patterns for the backward region. As noted, the former is just 
another example of the success of a Kirchhoff- type approximation 
in predicting the main beam of a radiation pattern,, whereas the 
failure to predict the patterns in the rearward region is largely 
the result of errors in the reflection coefficient resulting from 
these approximations. 

We conclude the chapter with some remarks concerning the 
technological relevance of the research described in this report. 

It has already been noted that the results of this study are of 
importance in the analysis of slot arrays radiating in the presence 
of a dielectric cover. As a matter of fact, it was pointed out in 
the first chapter that this motivated our research; in particular, 
we wished to assess the effect of the termination of the cover on 
the performance of the antenna. 

A dielectric cover is generally used for the purpose of imped- 
ance matching of the array (since the presence of the dielectric intro 
duces an additional parameter into the problem thereby providing a 
means to improve the scanning capability of the array by adjusting 
and optimizing the parameters), and/or as a protective cover for 
the array, i.e., as a flush mounted ra dome. The dielectric cover 
is of course intended to have a minimal effect on the pattern of 
the array. However, the work of Wu, Galindo, Lechtreck, Allen[18, 
19,20] and others have shown that if a surface wave is strongly 
excited in the cover, the energy radiated by the array can drop 
sharply thereby modifying the radiation pattern substantially. 

This condition is referred to as the forced surface wave resonance. 
Also, the use of dielectric cover for impedance matching may increase 
the frequency sensitivity of the array. Furthermore when the surface 
wave reaches the termination of the dielectric cover it radiates 
contributing to the pattern in a manner not intended in the original 
design. Using the diffraction coefficient derived here, a calcu- 
lation of this effect can be made. If the surface wave is strongly 
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reflected from the terminations of the dielectric cover (see Fig. 1), 
a large standing wave is created in the cover; this can have an 
undesirable effect on the input impedance and bandwidth of the array. 
This standing wave can be found from the surface wave reflection 
coefficient given here. 

The results of this study are also important to the design of 
surface wave antennas; a good review of the subject is given in 
Reference [21]. It is pointed out there that a surface wave guided 
along a uniform structure can radiate only at points of discontinuity, 
so that the terminal and feed radiation patterns are of interest. 

The configuration shown in Fig. 2 can be regarded as the termination 
of a flush mounted surface wave antenna, and the surface wave 
reflection and diffraction coefficients (which characterize the 
discontinuity in Fig. 2) may be used to advantage for analyzing 
the effect of the termination. It is seen from the results in 
this chapter that the magnitude of the surface wave reflection 
coefficient can be made small by decreasing the thickness of the 
dielectric panel; this in turn increases the level of the radiation 
field in. the forward region. In the case of a thicker dielectric 
cover,. it may be desirable to gradually taper the dielectric slab 
near its termination to reduce the magnitude of the surface wave 
reflected by the termination; this in turn improves the forward 
radiation characteristics of the end- fire surface wave antenna. 
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APPENDIX I 

, A WIENER-HOPF FACTORIZATION PROCEDURE j . . "'j: 

' • . ..u; 

I. In this section, a method is discussed for factorizing the' 
function L(s) into L+(s) and L_(s) which are respectively analytic 1 
in the U.H.P. and the L.H.P. As before, s = a + ix and L.H.P. - il, r 
refers to the region in the complex s plane for which x !i < T _ 

(where x_ is taken to be positive for convenience), whereas, U.H;P;' 
corresponds to the region in the s-plane for which x > x4(x+ is 
taken to be negative for convenience). L(s) has the properties 

. • .. tv;*' 

(A-T) L (s ) = L ( - s ) ' *- \ - ^ -I 


(A-2) L(s) is analytic and zonzero in the strip x + <' x < x_ 

(A-3) L(s) -v C A s v , when s -x» in the strip x + < x < x. 

C. and v are real constants. 


(A-4) Outside the strip of analyticity, L(s) may have simple 
poles, simple zeros, and branch points. 

As a consequence of (A-l), the poles, zeros and branch 
points occur in pairs, symmetrically with respect to 
the origin. For example in the case of the branch 
points s = tk (k = k-| + ik^; k^ << k-j and k^ > 0, k 2 > 0). 

The factorization of L(s) is achieved via a limiting operation 
on a related function K(s,a) which has only simple poles and simplb 
zeros; the factorization of K(s,a) is usually trivial. Specifically, 
K(s,a) is chosen such that 


(A-5) L(s) = 1 im K(s,a). 


K(s,a) need not be unique, even though its limit as a ->• is unique 
and equal to L(s). This technique of factorizing a function L(s) 
with branch points at s = +k through a limiting operation on a 
related function K(s,a) which has no branch points, was first de- 
veloped by Bates and Mittra [2]; they treated specific problems 
by this technique. Such a factorization technique generally leads 
to factors which are convenient for numerical calculation. More 
recently. Bates and Mittra[9,16] generalized their results in [2] 
to factorize a class of functions denoted here by L 0 (s) which 
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asymptotically behave as Cfts v e~ s for large s; C/\, v and b are 
real constants. L 0 (s) also satisfies the conditions indicated in 
(A-l), (A-2) and (A-4) for L(s). An integral representation for 
the Wiener-Hopf factors of L 0 (s) is initially conjectured by the 
authors in [16]; their conjecture is based on the form of the 
factors obtained earlier in [2]. In [16] the authors demonstrate 
their integral representation to be consistent with all the re- 
quirements on their factors, which lends credance to their fac- 
torization method. In this appendix, a derivation for the Wiener- 
Hopf factors of the function L(s) (with an asymptotic behavior 
indicated by (A-3)) is presented without requiring that the form 
of the factors be assumed as an initial step. The present develop- 
ment unlike that in [9,16] requires the introduction of a related 
factorization function K(s,a) as defined in (A-5); however, the 
final result for the factors is independent of K(s,a) even though 
the development of this factorization formula is based on initially 
introducing K(s,a) and later carrying out the limiting operation 
on, the factors of K(s,a) as a -*■ °°. In addition to (A-5), K(s,a) 
must be chosen such that it satisfies the following conditions; 


(A— 6 ) 

K(s) = K (— s ) 

t 

(A— 7 ) 

K(s) is analytic and non-zero in t+ < t < t_; 
elsewhere it has simple zeros and simple poles. 

(A— 8) 

K(s) -v C B s\ 

where C g is a constant. 

(A-9) 

lim K(s,a) -* 

a-x*> 

L(s), regularly in s 

(A— 1 0 ) 

lim K' (s,a) 
a-x» 

L' (s), regularly in s 


where the primes denote differentiation 
with respect to s. 


One may factorize K(s) into K+(s) and K«(s) where the + 
subscript implies analyticity in the U.H.P. and the - subscript 
implies analyticity in the L.H.P. as before. K (s) are given 
as [22] 
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oo 


s_ 

Pn 


-x (s) * o + r ) e rm 

(A-ll) K + (s) = [K(0)] 1/2 e 0 ^ ^ 


s 

oo e 9 

* (1 + — ) e 

n=l % 


n 


where K(s) has simple poles at s = ±q n and simple zeros at s = 1 p m . 
(Here, n = 1,2, 3 ••• and m = 1, 2, 3 •••). The function K_(s) is 
obtained from (A-ll) via 


(A-12) K_(s) = K + (-s). 


X 0 (s) is at present an undetermined entire function; it is chosen 
such that (A-8) is true. 


(A-l 3) X Q (s) = - X Q (-s); 

moreover X 0 (s) depends on 'a'; i.e, X Q = X Q (s,a). 

The representation in (A-ll) is converted to a contour integral 
as : 


( A- 14) K + (s) = 


[K(0)] 1/2 e ‘ X ° (S) exp --j^ 


J C UHP 


MW H 


K ‘ W d* 


The contour C UHp is indicated in Fig. A-l. 



Fig. A-l —The contour C UHp in the complex s plane. 
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The contribution to C UHp 


from C vanishes so that 

co 


(A— 1 5) 

1/2 _ ^ 0 (s) 

K + (s) = [K(0)] 1/2 e 0 
Define 



(A-16) lim X (s,a) = X(s) , 

a-*°° 


The limit as a-*» of ( A- 1 5 ) therefore leads to 


(fl-17) L + (S) = [L(0)] ,/2 e' X(s) exp 


& c 


■■ C * t) • tl *0 


Let L(s) have simple zeros at s = : (m = 1 , 2,3 •••) and simpli 

poles at s = tr p (n = 1 , 2 , 3 *** ). In addition L(s) has branch 
points at s = +k arisi ng from the presence of Vk^-s^ type terms. 
The proper branch for >/k2-s2 is the one for which Imvs^-k2 > 0: 


(the above is consistent with the development in Chapters II and III). 
One may now close the contour in (A-17) by Co, in the U.H.P. (with a 
deformation around the branch cut) and use the Cauchy Residue 
Theorem to obtain 
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The contour Cg r is shown in Fig. A-2. The contour -Cgr may be 
mapped into the contour Cw parallel to the real axis via. the 
usual transformation ^ = (k2-w2)l/2. Without going through the 
details of this transformation, the final result for the integral 
in the complex w-plane is given below: 


(A— 1 9) 


L + (s) = {L(0)} 1/2 




s 

I f(/? 

iVLi 

CM 

1 

CM 

J L(|k 2 - 

« 2 > 
s_ 

?? (1 

+ 

% 

m=l V 

x m/ 


CO f 

_ \ 

-s/ 

- (’ 

1 + — ] 

e 

w 




n 


The prime on L' (\/k 2 -w 2 ) now denotes differentiation w.r.t. w. The 

contour C, is indicated in Fig. A-3. 
w 3 
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Fig. A-3-- The contour C w in the complex w-plane. 


The contour may now be deformed to the real w axis. If 
has a zero at w = 0, the contribution from w = 0, which is 
by (1/2) v 0 must be included. 


(A-20) 


, . L' (Jk 2 -v?) 

lim w — ^ — — - L 

^ 0 L(jk 2 V) 


Hence (A-19) becomes 

(A— 21 ) L + (s) = [L(0)] 1/2 

s 




L(/k 2 -w 2 ) 


L(Jk 2 -w 2 ) 

denoted 
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Th e P.V. indicates that 1/2 (residues) of any other zeros of 

L(Jk 2 -w 2 ) encountered along the real axis must be included. The 
infinite integral in (A-21 ) may be re-expressed as a semi-infinite 
integral : 


V 


(A— 22 ) 

L + (s) = [L(0)] ,/2 e" X<S) " W 




m 


; ■ * W 


• f CO 

eXP ^ _ FriT P#V * [ 1n P + 

•'O' L ' \ 




{cf(w)+<jf(-w)}dwj 


where 


(A— 23) £(w) 


_ dw 


L(jkQ 

L Jk 2 V 


X(s) in (A-22 ) must be chosen such that L+(s) has the proper asymptotic 
behavior as indicated in (A-3). For specific problems, X (s ) may be 
determined from a knowledge of the asymptotic behavior of the infinite 
products and of the integral 


1 

2th 


P.V. 





w) + jf(-w)> dw 


which occur in (A-22); a typical procedure for the asymptotic 
evaluation of these terms as s-*» in the strip x+ < x < x_ is in- 
dicated in Sections II and III of this appendix. (Also see [2].) 

It is noted that the final form of the factors L+(s) indicated 
in (A-22) is not dependent on the choice of an appropriate closed / 
region function K(s,a); one may therefore proceed directly to the 
expression in (A-22) for L ± (s) without ever specifying a closed 
region function K(s,a). One notes however, that K(s,a) was 
introduced in order to arrive at the final expression for L + (s). 

In some cases it may be difficult to evaluate X (s ) in the ex- 
pression for L + (s) directly from the asymptotic behavior of the 
terms in (A-22)"; when this happens to be the case, it may be 
worthwhile to specify an appropriate closed region function K(s,a) 
for the problem in order that one may now obtain X(s) from the 

operation X(s) = lim X„(s,a). The function X (s,a) appearing in (A-ll) 
a-x» o o 

can usually be determined quite readily. 
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II. Application of (A-22) to the Factorization 
of il?7(Jk 2 -s 2 - ict)" 1 
The function L (s ) to be factored is 


(A-24) L(s) 



where a is a real constant. This function is encountered in the 
solution of the problem illustrated in Fig. 3. If one places a 
conducting boundary at x = a in Fig. 3, the function to be 
factorized becomes 


(A-25) K(s,a) = 


J k 2 -s 2 sin 


k 2 -s 2 a 


P ~ 2 


k^-s*" sin N k 2 -s 2 a +a cosJk 2 -s 2 


Clearly lim K(s,a) = L(s) which is also seen from physical con- 
a-x» 

si derations (i.e. moving the conducting boundary at x=a to ® yields 
back the original problem of Fig. 3). K(s,a) has no branch points 
in the s-plane. Both L(s) and K(s,a) satisfy the requirements 
indicated in Section I. Without re-deriving all the details in going 
from K+(s,a) to L+(s), one proceeds directly with (A-22). «J?(w) for 
this problem is 


(A-26) £M 


k 

L(,/k 2 -w 2 ) 


-iq 

w(w-i a) 


Also 


(A-27) v = lim w«£(w) = 1 
0 w->0 


and 


(A-28) jaw) + JSC-w) = - y 

VI + a 
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+ 2 2 2 

L(s) has no zeros; it has poles at s = _g (where k -g =- a ) and 
branch points at s = + k. Hence we may write L (s) for this problem 
via (A-21 ) as 


(A— 29 ) 






■oo 

•-00 


ln^l + 



s -ia dw 


Also, 


(A- 30) L_(s) = L + (-s) 


One knows from (A-24) that L(s) ^ CqS 0 , where Cp is a constant. 
Hence, L+(s) -v C J/ 2 ; one may now use this condition to find X(s) 
in (A— 297. 


The P.V. 


■oo 

in (A-29) may be written as 

J-co 


(A- 31) 
P.V. 




j (Residue at w=0) - 



The contours Cs, I and Cg r are shown in Fig. A-4. The integral 
over z gives a vanishing contribution. Only the 

In (7 + s/Jk 2 -w 2 ) 

term in the integrand can exhibit an algebraic behavior in s for 
L+(s) as s-*»; hence, one considers this term first. 


(A-32K1P.V. 

Zni 



It can be shown that 



is bounded and 


is a constant. 


Hence the 
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Fig. A-4 — The contours C s , £ and Cg r in the complex w plane. 


7 (' * k) 


dominates over I as s -* °° so that 
C n — 


\-33) exp^P.V. | ln(l + O ./(w)dw * 0(s 1/2 ) 


Utilizing (A-33) in the asymptotic behavior of L+(s) in (A-29), one 
obtai ns ' ; 


(A- 34) t 


L + (s) * [L(0)] V2 e' x(s) ' s/2k exp^ T RV.J 


- . k 2 -„ 2 


J?(w) dw 


/ c 1/2 

(since the ( 1 + r- ) factor in (A-29) behaves as s whereas 

/ C V-1 V K 

^1 + behaves as s ). 

Clearly the proper choice of X(s) such that L + (s) ^ constant as 
s-x», is 


,oo 

(a-35) x(s) --If -la p.v. 

J — C 


k 2 -„ 2 
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Incorporating the above result in (A-29) and converting the infinite 
integral representation of (A-29) to a semi-infinite integral 
representation given in (A-22) gives the final form for L + (s) as 


(A- 36) L + (s) = 





2 2 
w-wr 


dw 


L_(s) may be deduced from (A- 36) via (A- 30); alternatively, I (s ) 

may be obtained from the relation 


(A— 37) L_(s)=L(s)/L + (s). 

The numerical evaluation of the integral for L+(s) is carried out 
using a ten point Gaussian quadrature. 

III. Application to the Factorization of G(s)=(YY-|Sinh ^b)* 

• (Y-jSinh Y-jb+e^Y cosh Y-|b)"^ 


The factorization of G(s) = (yyt sinh Y-|b)(Y-|Sinh Yib + 
e^Y cosh Y-jb)"l is essential to the soluti on of t he problem 
discussed in Chapter III. As before, y = Vs^-k^, and Y] =J s 2 -k£. 

In this section, the factorization of G(s) will be carried out by 
proceeding directly to the expression for L+(s) in (A-22) without 
specifying an appropriate closed region function K(s,a). Here, 

G(s) corresponds to L(s) of Section I, and'9(w) will be defined 
to correspond toX(w) of Section I. Since the details for the 
factorization of G(s) are lengthy, only the important results leading 
to G + (s) will be given. Let 

(A-38) w’ =/ -k 2 +k^ + w 2 , where k and k^ are given in Chapter III. 
One May write^(w) as 


(A-39) 3(w) = ^ — - ~ W - = ^ In G(Jk 2 -w 2 ) 

G (J7V, 

where 


(A— 40). S(JkW) = iww ' sinw ' b 

-w* sin w'b-ie w cos w'b 
r 
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Hence, ^/(w) is given by 

4 2 3 2 3 

n w' sin w'b+ ie bw w' +ie w w' sin w'b cos w'b 

(A— 41 ) 3(w) 3—j E 5 

w[w' sin w'b + ie r w w' sin w'b cos w'b] 

Also 


(A— 42 ) » + 3( 


It is easily verified that v =1, since 

o 


(A-43) v = lim w^J(w) = 1 
w+0 


One may write G + (s) in terms of the above results via (A-22) as 



where the zeros of G(s) occur at s = m = 0, 1, 2 ••• and 3 0 = k d 
(see Chapter III for definition of B m ). The poles of G(s) occur at 
s = is (corresponding to the TM 0 surface wave poles). The choice of 
X (s ) is such that G+(s) behaves asymptotically as s*/2 when s-x»; the 
latter is a consequence of the fact that G(s) 'v 0(s) as s-x» in the 
strip within which G(s) is analytic. The asymptotic behavior of the 
integral 
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P.V. 


In 


0 + T J=Vp(w) + 3(-w) r 

V J k 2 -* 2 ' l J 


dw 


(or of it's alternative respresentation which involves the infinite 
integral 


P.V. 


r°° 

— 00 


In 





occuring in the expression for G+(s) in (A-44) is ascertained in a - 
manner identical to that in Section II of this appendix (see (A-32) 
and ( A- 3 3 ) ) . The infinite product. 



occuring in (A-44) is replaced by 



so that it's asymptotic behavior is readily obtained via [9]; thus 


(A- 


ist) 

IT 


■45) ~ (l + §-) e n 

n=l V 


sb isb 
.s-^e 2 e ' 


r o* 1 




arg t t w. 


One notes that changing the exponential factors in the infinite 
product is justified since the negative of these exponential factors 
must also be present in a similar representation for G_(s); these 
exponential factors are required to ensure the convergence of the 
infinite products. One notes that ro in (A-45) is defined as 

Euler's constant = 0.5772156649. 

It is easily seen that the proper choice for X(s) which ensures an 
algebraic growth for G+(s) as s-k° is 


(A-46) X(s) =M[ ro + ln (|k)] 


sb 


+ I(s) - 5 
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where 


(A— 47) 


I(s) 


oo 

dw 

b 2f(w) +^(-w) 

Jk 2 -„ 2 

2-rri 


Hence G + (s) is 


(A-48) 


) 1/2 (s+k d ) 

W f- 

f g(s,w) - £ + G(w) dw 
. e 0 L J e P(s) 



G + (s) = g] /2 (0) 


(s+k 


(’ 


where 


(A— 49) 


g] /2 (0M 


i sin k d b 


11/2 


-k d sin k d b - ie k k d cos k d b 


I 

f 


(A— 50) g(s,w) = In (l + — ^ V 

^ I 0 o j 


yy j k w 


(A-51) G(w) = 


. , e w 2 b + f — r — w') g sin w'b cos w'b 

b^ _r ^w J r . 

71 77 w' 2 sin 2 w'b + g 2 w 2 cos 2 w'b 


and 


- r - I(s) • 

Finally, G_(s) is identical to G+(-s). The representation in (A-48) 
is not suitable for numerical computations, because I (s ) in (A-52) 
is not well behaved. If one consideres the contribution to I(s) 
over the range w 0 < w < » where w 0 is an arbitrarily large number, 
then one may accurately approximate the integrand as follows: 


(A-52) P(s) = 


isb 


+ In ^ 



2tt 
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J Tj- 


b 1t M +*c((w) 

“ A • 


z^r 


- T- -(' 

1W IT \ 


l+( e 2 -l )C0S 2 Wb 


) 


The integration of the above from w 0 to « would converge if the term 

(l e r 

^ l + (e 2 -l )cos 2 wb 

was purely oscillatory (i.e., behaved as sin w); however, the 



1+(e 2 -1)cos 2 wb 

\ 

term in the brackets oscillates about a non-zero mean value equal 
to 


e r + 1 

2e 


(which will equal unity only if e r =l), thereby implying that the 
integral I (s ) of (A-47) is non-convergent for There are no 

other terms in G+(s) which can compensate for the divergent 
behavior of I(s), thereby making the above representation for 
G+(s) unsuitable for numerical computations. This factorization 
procedure has been applied to a number of examples in the past 
without difficulty; hence, the present result is unexpected. In 
this sense it is a peculiar case. In Appendix II we find G + (s) by 
a formal factorization procedure. 


104 



APPENDIX II 

EVALUATION OF THE WIENER-HOPF FACTORS OF G(s) VIA THE 
FORMAL FACTORIZATION PROCEDURE 


I; It was indicated in Section III of Appendix I that the limiting 
type factorization procedure of Appendix I although ordinarily pre- 
ferred over a formal factorization procedure for convenience in cal- 
culating the factors, strangely leads to a representation for the 
factors of G(s) which fail to converge. Therefore, the factorization 
of G(s) is carried out via the formal factorization scheme described 
in this appendix. Since the well known texts on the subject of 
Wiener-Hopf methods such as those by Noble[19] and Mittra and Lee[9] 
deal with the development of formal factorization formulae, only a 
summary of the results will be given and then applied to the factori- 
zation of G(s). 

Consider the complex function F(s) of the complex variable s with 
the following requirements: 


(A-53) F(s) is analytic, regular and non-zero in the strip 
x+ < t < t_ (where t * Im s; s = a + it). Figure 10 
indicates the strip of analyticity for F(s ) . 


(A-54) F(s) -> 1 uniformly . as |o| » in the strip 

T + .< T < T_. 

It can be shown that F(s) may be factorized in the folliwing 
manner: 


( A- 55) 
where 

( A- 56) 


and 
(A— 57) 


F(s) = F + (s) F_(s) 


F + (s) = exp 


1 

27 ri 


oo+l T. 


-oc+i i _ 




T >Im S>x >T, 
— a « 


F_(s) = exp 


1 

2TTi 


oc+i x 


■oo+i X, 


dz 


In F(z) 


z-s 


x + <Im S<T^<T_ 
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The contours of integration for the integrals in (A-56) and (A-57) are 
shown in Fig. A-5(a). F + (s) is regular and non-zero in the upper half 

s plane defined by t > x+, whereas F_(s) is regular and non-zero in 
the lower half s-plane defined by t < x_. It is possible to evaluate 
F+(s) for Im s < x a and F_(s) for Im s > xb by deforming the contours 
of integration in Fig. A-5(a) to enclose the point s as indicated in 
Fig. A- 5 ( b ) (for s = s-| with Im s-| > x b and for s = s 2 with Im s 2 <x a ). 

G(s) of (94) may be factorized via the formulae of (A-56) and 
(A-57). G(s) is re-written below for convenience. 

(A-58) G(s) = (yy-j sinh Yi b )(y-| sinh Y-] b + e r Y cosh Y-| b ) _1 
where 

(A-59a,b) Y = v/s 2 -k 2 = -ijk 2 -s 2 , Y] = J s 2 -k 2 , and 

k and kj appearing in (A-59a,b) correspond to the wave numbers in 
free space and in the dielectric slab (of thickness b), respectively. 
The zeroes of Yisin hY]b in (A-58) occur at s = tk^ and at s = ±i|Bnl 
(where n = 1,2,3,***)« Since the problem is restricted so that only 
the TM 0 surface wave propagates, the only zeroes of (yi sinh Y ]b + 
e rY cosh Yl b ) in (A-58) are at s = ±g. The branch points of G(s) occur 
at s = ik and at infinity. The branch points, zeroes, and poles of 
G(s) are indicated in Fig. A-6. For convenience of analysis, k and 
kd are made complex by allowing them to possess small imaginary 
values (this corresponds to introducing a slight loss in the medium 
as indicated in Chapter III) so that 

k = k-j + ik^ f with k 2 « k-j \ and k d = k 3 + ik^ 

Vand k-| , k 2 > OJ 

this allows G(s) to be analytic and non-zero within the strip, defined 
by x+ < x < x_, and |a| < ”, in the complex s plane. x+=-min(k 2 ,k 4 ) 
and x_=min(k 2 ,k 4 ) ; T± -* 0 as k 2 and k 4 both approach zero. Figure 
A-7 indicates the locations of the poles, zeroes and branch points of 
G(s) when k and kd are allowed to be complex. Since G(s) ^ 0(s) as 
J s | -* oo in the strip, it is not possible to factorize G(s) directly 
via (A-56) and (A-57) since the condition in (A-54) is violated. It 
is therefore necessary to re-express G(s) as 


with k^ << k 3 \; 
and k 3 ,k 4 > 0 J 


(A-60) G(s ) = G-, (s ) G 2 (s) G 3 (s) 
where 

(A-61) G ] (s) eJ s 2 -k 2 = -ijk 2 -s 2 
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Fig. A-5 — Contours of integration for the Wiener-Hopf factors. 
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k 0 


a? Re s 


Fig. A- 6 - 




-Poles, zeros and branch points of G(s) in the 
complex s plane. 

Ims (OR Imz ) 


CONTOUR OF INTEGRATION 
Gj (s) IN THE 
Z - PLANE \ 



u 1 o i,*. * + T \ I 

~ k d P K \-0i» \ POINT Z = S T q 

\ \ IN THE 

\P Z » \ COMPLEX 

\ \ Z - PLANE 

V« CONTOUR OF INTEGRATION 

1 I FOR G+(s) IN THE 

I ' Z -PLANE 

* s = <r +ir PLANE 

OR 

Z = Z, iz 2 PLANE 

Fig. A-7 —Zeros and singularities of G(s) when k and kd are complex, 
and the contours of integration for G^(s) in the z plane. 




(A— 62 ) 

and 

(A-63) 

with 

(A-64) 


g 2 (s) 5 (s 2 -k|)/(s 2 -e 2 ) 


sinh„ 

Ms) s -= 

J i o 


s2 - k d b s 2 - e 2 


sJS K, 


g(s) 


g(s) = >/s 2 -k^ sinhjs^-k^b + eJT-7 coshj s 2 -k ^b 


G](s) and G 2 (s ) may be factorized by inspection, or via the 
procedure given in Section I of Appendix I. Thus, 

(A-64a,b) G 1 (s) = g|(s) G"(s); g{(s) = (-ik) 1/2 (l + f-) V2 = ^(-s) 


and 


(A-65a,b) G 2 (s) = G+(s) G~(s); g£(s) = 


0 J 2 

('•H 


= G'(-s) 


G 3 (s) is analytic and non-zero in the strip (t+<t<t_; |o|->«>); 
furthermore, G 3 (s) -> 1 as I s | -* <=° in the strip. Thus, G~(s) may be 
factorized into G^Cs ) G^(s ) via (A-56) and (A-57). 

(A— 66a ) 

G^s) = exp k + (s) = exp wU 


00+1 T a In G~(z) ' 

J dz 


- 00 +i X. 


z - s 


T_>Im S>T a >T + 


(A-66b) 


G’(s) 5 exp k_(s) = ex P | " 2rf J_^ u 


oorf i x 


b In G 2 (z) 


z - s 


dz 


|x + <Im S<T|^<T_ 
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The quantities k+(s) and k_(s) are the terms which appear inside the 
brackets on the R.H.S. of (A-66a) and (A-66b), respectively. The 
contours of integration and the location of z = s in the complex z 
plane for (A-66a,b) are indicated in Fig. A-7. One may express 
k+(s) as 


(A— 67) 


M s ) = 2 ui 


r 



JO+i t 

a 


In G 3 (z) 
z - s 



In G 3 (z) 
z + s 



One now lets k£ and k 4 approach zero (reducing the loss in the medium 
to zero); this allows one to express k + (s) above as 


(A-68) 


k + (s) 



In G 3 (z) 



dz 


where the contour of integration C+ is indicated in Fig. A-8(a) 
whenever s is real and positive. A similar representation for 
k_(s) when k2 and k 4 are set to zero (i.e., when the strip of 
analyticity for G 3 (s) shrinks to the real z-axis, except for the 
singularities at z = k and z = s) is given by 


(A— 69 ) 



In G 3 (z) 



dz 


where the contour C_ is depicted in Fig. A-8(b) for positive real s. 



V 


s 

■M? * 

CONTOUR C. 


TO +CO 

i 


(a) 


A 


K ' s 

^CONTOUR C_ 


(b) 


Fig. A-8 — Contours C + and C_ associated with G 3 (s) and G^(s ) for real s 
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It can be shown that the contribution to the integrals in (A-68) 
and (A-69) arising from the branch point at z = k is vanishingly 
small; hence, these integrals may be expressed as principal value 
(P.V.) integrals over a semi-infinite range along the real z axis 
after taking into consideration the effect of the pole at z = s 
(by including one half of the residue arising from this pole). 
Therefore, 

C» In G~(z) 1 

(A- 70) k + (s) = f- P.V. — T dz + f In G,(s); for s>0 

Jo z - s 6 

and 

In G.(z) , 

(A-71 ) k (s) = P.V. — 9 — ■ — - 9 dz + y In G,(s); for s>0 

11T JO z - s 

Thus, one may factorize G(s) into G + (s) G_(s) where 

(A- 72a) G + (s) = g|(s) G^Cs) G^s) 

and 

(A- 72b) G_(s) = G“(s) G“(s) G“(s) 

with Gf(s) and G^(s) indicated in (A-64b) and (A-65b). G|(s) in 

(A-72a,b) are defined in (A-66a,b) with k+(s) and k_(s) given by 
(A- 70) and (A-71), respectively. Equations (A-70) and (A-71) are 
valid for s >0 ; when s<0 (i.e., s=-| sp one may employ the relations 


(A-73a,b) k + ( - |s | ) = k_( |s | );k_(- |s () = k + ( |s 

The relations indicated in (A-73a,b) follow from the fact that 
G 3 (- |s I) = G 3 ( \s |). Clearly, 


(A-74a,b) Gjdsl) = GjHsI), and - |s | ) = G'( |s | ) , 

whenever (A-73a,b) are valid; therefore, it follows via (A-72a,b) 
that 

(A- 75) G + ( |s | ) = G_(- |s | ). 
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The factors of G^(s) when s is complex may be evaluated directly via 
the, formulae in (A-58) and (A-57) with the corresponding contours of 
integration indicated in Fig. A-5. The relationships 

(A- 76) G^Cs) = G~(-s) '..f 

and 

(A-77) G + (s) = G_(-s) 

are valid even for complex values of s. 

Specifically, one needs to evaluate G + (s) at s = B,,s.= kn, . • 
s = 0 n (='i|3 n |, n = 1,2,3 • ••) and s = k sin <p (with j 7t/2), i 
in order to construct the surface wave reflection and diffraction." 
coefficients pertaining to the configuration of Fig. 2; however, 
before proceeding to these cases, it is desirable to re-write 
(A-70) and ( A- 71) in a form suitable for the numerical evaluation 
of the integrals. 


The quantity 


t j- P.V. 

1 TT 


In G 3 (z) 

~2 2~ 
z - s 


dz 


in (A-70) and (A-71) may be re-expressed for the purposes of numerical 
computations as ' 


(A- 78) 


/ 


+ 4- P-V. 

- 1 IT 


+ s_ 

iir 


= S 


+ s 


Tim 

e-K) 


linf< 


k-s 


1ir Uo HO 


S -5 


fS -6 j 

S-E 

rS +6 

+ 

•'k+E J 

+ 

S -6 

S+E 

■S-E 

rS +6 , 

k-E 

+ 

U-S J 

s + e + J 

S +6 


S +6 


,for s>k>0 




for 0<s<k 


One notes that the integrand in (A-78) is real whenever z>k and s is 
real. It is desirable to split the integral as shown above since it 
allows one to substantially reduce any numerical roundoff errors 
arising from the numerical evaluation of the integral in the neigh- 
borhood of the pole at z = s; the integral in the deleted neigh- 
borhood of the pole is then evaluated very accurately in closed 
form as discussed below. 
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G 3 (z) is analytic at z = s as well as for all positive 
real values of z excluding the point z = k; hence. In G 3 U) in the 
integrand of (A-78) may be accurately approximated in the irrmediate 
neighborhood of z = s by it's two term Taylor Expansion as 


(A- 79) In G 3 (z) £ In G 3 (s) + 


oiyi? G 3 (z) 


(z - s) 


-j z = s 


G'(s) 

" ln G 3^ s ^ + GJTsT 


(z-s). 


The integral from z = s-s to z = s + 6 excluding the point z 
which occurs in (A-78) may be written via ( A- 79 ) as 


(A-80) ± ~ lim 

llT e+0 


■s-e In G 3 (z) 


L's-6 


z 2 - s 2 


dz + 


•s +6 In G.-(z) 
dz J 


1 S+E 


2 2 
z - sr J 


= + f— lim ln G-(s)' 

- ” ^0 3 


S-e 


dz 


J S-6 z 2 - S 2 


■S+6 


dz 


Go(s) f fS-t 




dz 


S-6 


z + s 


rS+6 

U+z 


S+e z 
'v 


dz 


z + s 


The extent of the interval 26 about z = s is chosen such that the 
R.H.S. of (A-79) is a sufficiently accurate representation for 
In G 3 (z) within the interval. After some straightforward manipu- 
lations (A-80) becomes 


(A- 81 ) 


t ~ lim 
lTr e+0 


S-e 


In G,(z) 


dz + 


= i lim 
l7r e+0 


= + ]_ 
TiT 


L: S-6 Z - S 

In G 3 (s) 


s +6 In G 3 (z) 


2 2 
s+e z - s 


dz 


2 s‘ 


(6-e) + 


G ^(S) 

sG 3 (s) 


(6-e) 


In G 3 (s) GJ(s) 


2 s 


+ W 


Incorporating (A-81) in (A-78) yields 
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GO 

(A- 82) + 4-P.V. ln G 3 (z) 

\ TT / /•> o O 


10 z 2 -s 2 


dz = 



fk-e 

rS-6 

+ 

+ 

Jo j 

k+e 


fS-6 

■k-e 

+ 

+ 

Jo J 

S+6 J 


S+6 


k+e 


+ T_ 

1 IT 

s>k 

+ 1 _ 

i 7T 

J s<k 
s> 0 


ln G 3 (s) GJ(s) 


2s 

> 3 ( 


+ w 


ln G,(s) G'(s) 

+ b 


2 s 


G 3 (s) 


The result given in (A-82) may be directly employed in the integrals 
on the R.H.S. of (A-70) and (A-71) to evaluate K+(s) and K_(s), 
respectively when s is real. When s has a non-zero imaginary part 
(s is complex), the representation for k+(s) and k_(s) is taken 
directly from (A-56) and (A-57) respectively, together with the 
corresponding contours of integration indicated in Fig. a-5. 

II. Evaluation of G^ ( 3) : The factors of G 3 for s=B where B>k>0 are 

obtained via (A-70) and (A-82); in particular,. G^C B) equals e^ + (^ 
where 


(A-83) k + (B) = 


1 TT 


lim|f k £ 1n G 3 (z) dz + f 


U° z 2 - b 2 


B-6 In G 3 (z) 
•'k+e z 2 - B 2 


dz + 


j* 00 In G 3 (z) 
_ _ 

3+6 Z - 3 


. I In G,(B) G'(b) 

dz + ^ + 7 6 + 


T IT 


2 B 


Go(b) 


1 


+ J In G 3 (b) 


The term, ln 63 (B) appearing in (A-83) is calculated via 

/ 2 bV? > 


(A- 84) ln G-(B) = In 


lim G~(z) 
L z+3 - 


In 


-(^- k d)J^-k Z > 


where ^ is the residue defined by lim (z-b)G(z) which was introduced 
in (113). The quantity 63 (B) in (A-83) is presented later in (A- 88 ). 
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One notes that the ratio (z -3 )/g(z) occuring in the ex- 
pression for G3 (z) (see (A-63)) is sensitive to any slight inaccuracy 
involved in the numerical evaluation of 3 whenever z is in the 
neighborhood of 3 . Thus, it is apparent that the numerator and the 
denominator of G3(z) may not simultaneously attain a zero value at 
z=3 if a sufficiently precise value of 3 is not available. In order 
to make the integral in (A-83) less sensitive to any small errors in 
3 , the integral from k+e to 3-6 and also from 3+6 to 00 are further 
split as follows 


(A- 85) 


(-3- 6 

* OO 

'3-T 

r3- 6 , 

r3+T 

+ 

= 

+ 

+ 

+ 

•’k+e 

3+6 

k+e 

3-T J 

3+6 ' 


and the term g(z) which appears in G3(z) for the integrals over 3-T 
to 3 - 6 , and over 3+6 to 3+T, is replaced by it's four- term Taylor 
expansion about z=3. The interval 2T about z=3 is therefore chosen 
such that g(z) is accurately represented by four terms in a Taylor 
expansion for g(z) about z=3. In the interval 3-T < z < 3+T, g(z) 
is therefore given by 


(A— 86 ) g(z) 4 g(e) + g'(8)(z-8) + j- g"(6)(z-e) 2 + jg"' (B)(z-B) 3 

- Cg'(s) +ig"(B)(z-8) + ^r(e)(z-B) 2 ](z-B) 

since g(l3) = 0. T he quan titie s g'(z) an d g "(z ) are given below in 
terms of y(=Jz 2 -k^ = - ij k 2 -z^) and rj (njz^-ks) which have been 
defined earlier; these quantities are evaluated at z=3 and employed 
in ( A— 86 ) above. 


(A-87a) ’ g(z) = y-j sinh y-jb + e r y cosh y-j b 


z_ 

Y 1 


(A-87b) g'(z) = 77- sinh y-|b + zb cosh y-jb + e r “ cosh Y-| b + 


+ e„ z b sinh y. b 
r y. '1 


(A-87c) g"(z) = 


sinh y.b 2 


V z‘ .... „ , z 2 b 


Y 1 


3 

Y 1 


sinh y-|b + cosh y-jb + 


z 2 b 2 


+ b cosh y-jb + sinh y-jb 4- 
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(A-87c) cosh y,b 2 ■ 2. 

(Cont). + e- — 1 — — — - e 7 ? cosh y -I b +‘ 2e„ - — sinh yvb + 

,,-r . jY , : r . 3 .. - 'I.; / . :■ T yyy • 'I 

•:-f •• \ ’Tv 2 2 ‘ ; ; '• '' * 

^ sirih yvb + e cosh 

■. ■ . ! ' r TJ - . J- y|. : ■ ■ ... 


*1 


b - 


- •£• Y“^ sirih y.b. 

< . -i ■ » . >1 • • ,0 *. - I . 

■ ' - Y i • 


The expression for g"'(z) is rather lengthy, but can be obtained by 
a straightforward differentiation of g" (z); hence, it will not be 
presented here. The quantity G-j(b) in (A-83) may be expressed in 
terms of g 1 ( b ) and g"( 3 ) as 


(A- 88 ) GJ(e) = _____ 

sin Jk 2 -B 2 b 2 B 2 b cos^K 2 rB 2 b ' 


g* (e) 




sinjk2-B 2 b 


L M 


- g"( 6 ) • b 


sin 


2~J 


(k 2 d - 6 2 ) 


(k 2 - B 2 ) 


1,2 J 

Jk d -B 


i4|- 

Jk 2 -B 2 


B b 


[g'(e )] 2 


III. Evaluation of G 3 (k d ): The evaluation of G 3 (k d ) where k d >B>k 

proceeds in a manner similar to that for GjfB) discussed earlier. 
Specifically, G 3 ( k ^ ) = exp(k + (k ( j)) where 

r k-e In G,(z) kj ' 6 

(A- 89) k + (k d )=^ 1 



z 2 -k 2 
z K d 


dz + 


"d~ u In G 3 (z) 


k+ e 


z 2 -k 2 

Z K d 


dz + 


f In G,(z) 

° A- 

r 

+ i_- 

In G 3 (k d ) ^ G ^(k d ) ^ 

2 2 az 
V 5 2 - k d, J 


1 7T 

[ 2k d ,.WJ 


rs + 


+ 7 T " G 3< k d> 
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and 


b(k 2 -B 2 ) Jb(k 2 -B 2 

^ W = — F= •• ln W = lnJ 


)] 


e 

e W K d K 


k 2 -k 2 


rg d 


Also 

(A-91 ) GJ(k d )‘ 


( 


< \2k 

:Jk 2 -k 2 . 


d b - 2k d b • G 3 (k d ) 


r "d 


k d b(k d" 02} 3 2 2 

- V - S - - k d b < k d-8 ) 

(k 2 - k 2 ) d d 


' 2 
k.b^ 


3 * G 3^ k d^‘ 


Since B lies in the interval k + e < b < k^ - 6, it is desirable to 
re-express the integral over k + e to k^ - 6 in (A-89) as 


(A-92) 


k d" 6 
■ k+e 


B-T 


b+t 


k+e J$-T 


K d -6 

B+T 


The g(z) term in G 3 (z) for 


•B+T 

• b-t 


In G 3 (z) 



dz 


of (A-92) is replaced by it's four term Taylor expansion about z = b as 
indicated earlier for G+(b) for the purposes of making the calculation 
of G 3 (z) insensitive to any small differences which may exist between 
the calculated value of b and the true value of B. 


IV. Evaluation of G 3 (i|$ n |): G 3 (i|e n |) = exp k + (i.|B n |)= G 3 (g n ) 


where 




as indicated in Chapter III, and 


117 



(A— 93) K + (i |3 n | ) = 




1 IT 


- 

r 

lim * 


_g-*0 

. 


k-e In 63 ( 2 ) 


z2+ I8 n | 2 


dz + 


f3-T In Gg(z) 
ktE z 2 +|S n | 2 


dz+ 


•p+T In G 3 (z) 
J P“T z 2 + ! 3 n 1 2 


dz + 


r ln g 3 (z) 

b+t ^7 


dz> 


The term g(z) in G 3 (z) which occurs in the integral over p-T < z < 3 + T 
in (A-93) is approximated by it's four term Taylor expansion about z = 3 
as done previously for G"§(e ) and G^kj), and for similar reasons. One 
notes that the integrand in (A-93) possesses pole singularities on 
the imaginary z-axis unlike that for k+( 3 ) and k+(kd) for which the 
singularities occur on the path of integration. 


V. Evaluation of G 3 (k sin $): Since |<j>| < tt/2 (see Chapter III for 

detail it follows that |k sin <j>| < k. Thus, G^(k sin <|>) = exp 
k+(k sin <|>) • = exp k+(±|k sin <j)|), for <f> £ 0. It was indicated in 
Section I of this appendix that one can relate k+(-|k sin <}>|) to 
k_(|k sin ((>]) by the equality 


(A-94) k + (-|k sin <j> | ) = k_(|k sin cj>|) . 


The quantities k + (|k sin <{> | ) are given by 
(A— 95 ) 


k + ( |k sin 4 >|) 


|k sin <p 

i TT 


1 im 

e^o Mo 


[k si n 4> (-6 In G 3 (z)dz 


2 2 
z -|k s i n 4 > | 


rk-e In Gg(z) dz rP-T ln G 3 (z) dz ' 

|k sin <p | + 5 z 2 - 1 k sin <}>| 2 ^k+G z 2 -|k sin cj >| 2 


3 +T In G 3 (z) dz 


2 ? 
8 -T z - |k sin <f>| 


dz + 


r°° In G 3 (z)dz 


2 2 
8 +T z - 1 k sin f| 


J - * 
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(A-95) -In G,(|k sin 4 > ] ) G ' ( | k sin <j>[) 

(Cont) + {-"S + — [6 + 

77 2 1 k sin <j> | G~( jk sin <p | ) 

v 6 J 

+ 1 In G 3 ( |k sin <j> j ) , for 4 > > 0 


(A-96) 


k_( |k sin <p | ) =' 


f r jk si n«f> | - 6 In G~(z)dz 

lim ‘i -? ? 

£-*0 I 'o z 41 - | k si n<j> | 


In G 3 (z) dz rB-T In G 3 (z) dz 


■'Ik sin <f> |+6 z 2 -|k sin <}>| 2 + ■'k+e z 2 - | k sin <f>| 2 


f B+T In G^(z) 

+ o dz + 

•'b-T z - |k sin <j»p 


« In G 3 (z) 

B+T z 2 - | k sin <j>| 2 


dz> - 


1 J" 1n G 3 ^l k sin ♦!) G 3 d k sin <H)|_ 

T7 j 2 |k sin *| + G ( | k sin + 


+ j In G 3 ( ) k sin <j> | ) * for </> < 0 . 


In order to evaluate G^jk sin <{> } ) one must set z = | k sin <p | in 
(A-97) r , o 


«3< z > ’ 


sinh Y 1 bl( Yl sinhy 1 b+ e r Y coshy-, b) (2z )- (z -B )fz) 


/sinhy,b Coshy,b 

•I + b cosh y,b +e + e 

\ Yi 1 r Y r 


yb sinhy-jb 


( Y] sinh y-jb + e r Y cosh y-, b)" 


119 



(A— 97 ) 
(Cont) 




Finally, the term g(z) in G 3 (z) which occurs in the integral over 
6-T to 8+T in (A-95) and (A-96) is approximated as in (A-86) for 
the purposes of making the numerical evaluation of the integral 
less sensitive to any small differences which may be present in 
the calculated value of 8 from it's true value. 
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